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Preface 


The  purpose  of  thia  atudy  waa  to  investigate  tha  af faotivanaaa  of 
uaing  an  alaotroatatio  precipitator  aa  a  maana  for  reducing  harmful 
lavala  of  radon  progany  in  tha  homa.  Thia  ia  an  important  atudy  baaauaa 
mora  and  mora  homaa  may  ba  axparianoing  buildup  of  dangacoua  lavala  of 
radon  and  radon  progany  dua  to  axtanaiva  afforta  to  aava  anargy.  By 
aaaling  up  tha  homaa,  air  axchanga  rataa  ara  drastically  raducad  thua 
oauaing  radon  daughter  oonoantrationa  to  build  up  in  tha  homa.  During 
thia  invaatigation  1  will  look  at  a  method  whioh  may  ba  effective  in 
raduoing  theaa  dangaroua  lavala. 

In  preparing  for  and.  carrying  out  tha  experimenta  in  thia  thaaia  X 
received  a  great  deal  of  help  and  guidanoa  from  ethara.  X  am  deeply 
indebted  to  my  thaaia  adviaor,  Dr.  George  John,  for  hia  guidanoa,  knowl¬ 
edge  and  ocoaaional  kick  whan  tha  need  praaantad  itself.  I  would  also 
like  to  thank  Bob  Handrloka  for  hia  support  and  aaaistanea  with  muoh  of 
tha  equipment.  A  apaoial  thanks  goes  to  Carl  Shortt  and  the  members  of 
tha  AFXT  fabrication  shop  for  their  asaistanoa  in  building  tha  Radon 
Test  Chamber.  Additionally,  X  would  like  to  thank  Dr.  Hagee  and  hia 
aasooiate,  Dr.  Jenkins,  of  Monsanto  Mound  Facility  for  asaistanoa  in 
calibrating  tha  large  Lucaa  cells  used  in  radon  gas  measurement.  Finally, 
I  want  to  thank  my  wifa  Ann  for  her  untiring  support  and  devotion  during 
thoaa  momenta  when  I  wondarad  if  it  waa  all  worth  tha  effort  and  pain. 


David  R.  Little 
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Aba tract 

The  effaotivenaaa  of  uaing  an  alaotroatatio  praeipitator  aa  a  maana 
for  reducing  harmful  lavala  of  radon  progany  In  tha  homa  waa  avaluatad. 

A  oommarcially  availabla  praeipitator,  manufactured  by  tha  Honeywell 
Corporation,  waa  uaaa  during  tha  eouraa  of  thia  atudy.  Tha  apaoifio 
modal  uaad  waa  tha  Honaywall  Elaotronio  Air  Claanar  modal  numbar  F50A1009. 

Daughtar  oonoantrationa  wara  maaaurad  by  tha  modifiad  Taivoglou 
mathod.  Samplaa  wara  oollaotad  on  a  2-inch  millipora  filtar  and  alpha 
amiaaiona  wara  maaaurad  with  a  ZnS(Ag)  acintillator .  A  aampla  aollaotion 
tima  of  3  minutaa  waa  uaad.  Sampla  counting  intervale  of  2-3,  6-20,  and 
21-30  minutaa  aftar  aampla  oolleotion  wara  uaad. 

During  thia  atudy  air  aamplaa  wara  mada  uaing  tha  blowar  fan  and 
filtara  alona  with  no  powar  to  tha  alaotroatatio  praoipitator,  and  with 
tha  alaotroatatio  praoipitator  anargisad.  Tha  raduotion  in  tha  working 
laval  aa  a  raault  of  uaing  tha  blowar  fana  and  filtara  only  waa  73  paroant 
With  tha  alaotroatatio  praoipitator  anargiaad  tha  raduotion  laval  roaa  to 
90  paroant.  It  ia  tharafora  eonoludad  that  tha  alaotroatatio  praoipitator 
ia  an  affective  maana  for  raducing  radon  progany  oonoantrationa  in  tha 
homa. 
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ANALYSIS  OF  RADON  AND  RADON  PROGENY  IN  RESIDENCES i 

FACTORS  THAT  AFFECT  THEIR  AMOUNTS  AND  METHODS  OF  REDUCTION 

I .  Introduction 

Purpose 

The  purpose  of  thia  thaaia  waa  to  study  and  quantify  the  ef-'ective- 
naaa  of  using  an  eleotrostatle  air  olaanar  in  reducing  tha  radon  progany 
ooncantrat ions  in  homaa. 

Background 

Inhalation  of  air  containing  radioactive  radon  gaa  or  progany  nay 
result  in  an  axpoaura  of  tha  raapiratory  traot  to  harmful  lavala  of 
ioniiing  radiation.  Prolonged  axpoaura  to  thia  radiation  has  baan 
attributed  to  a  high  inoidanaa  of  lung  oanoar  oaeaa  among  mine  workers 
in  tha  United  states  and  Eastern  European  oountriaa  (1iS6>7i98),  Evan 
though  knowledge  of  tha  potential  radiation  ha card  exists,  vary  few 
efforts  have  baan  direated  toward  tha  eontrol  and  raduotion  of  this 
hasard.  Efforts  to  aontrol  this  hasard  are  hampered  because  of  tha 
leak  of  understanding  of  the  complex  relationship  between  radon,  its 
progeny,  and  their  environment  (23il.l).  It  was  believed  originally 
that  the  lung  cancers  in  mine  workers  were  dus  to  high  levels  of  rad^n 
gas i  however,  W.  F.  Bale  demonstrated  that  under  ail  practical  circum- 
stances  the  dose  to  critical  tissue  in  tha  lung  area  was  due  predomi¬ 
nate^-  to  tha  short  lived  daughters  of  Ra-222  (Po-210,  Pb-214,  Bi-214, 
and  Fo-214)  which  are  present  in  radon  laden  air,  rather  than  to  the 


radon  gas  itself  ( 4 1  323 )  . 


Not  long  ago  It  was  belJwved  that  the  only  parsons  who  ware 
exposed  to  signif iaantly  high  levels  of  radon  and  its  daughtara  wara 
undarground  mina  workers.  Tha  eonoantration  of  radon  gas  was  expsctsd 
to  be  high  because  of  tha  larga  quantity  of  uranium  or  radium  oontant 
oraa,  tha  larga  surfaea  araa  for  amanation,  and  tha  raduoad  vantilation 
(23 i 1.1).  In  raoant  timaa  many  articlaa  hava  baan  publiahad  expressing 
growing  oonoarna  toward  tha  ganaral  public's  exposure  to  haiardous 
levels  of  ionising  radiation.  Initially  oonoarna  wars  diraotad  toward 
thoaa  haaarda  producad  by  man  (i.e.,  nuolaar  powar  plants,  uaa  of 
nuolaar  waapona,  and  madical  application*  of  radiation),  Mora  raoantly 
though,  natural  aouroas  of  radiation  hava  raoaivad  inoraasad  attention , 
Tha  moat  racant  artiolaa  hava  baan  oonearnad  with  tha  magnituda  and 
variation  of  indoor  radon  progany  inaida  raaldantial  houaaa  (24|27|32» 

33 1 38 ) .  Thoaa  raporta  hava  shown  that  Indoor  exposures  to  radon  progany 
may  ba  several  timaa  graatar  than  that  experienced  outdoors.  Thaaa 
alavatad  conoantrations  ara  usually  asaooiatad  with  homos  which  hava 
baan  oonatruotad  on  top  of  soil  with  high  uranium  or  radium  oontant, 
resulting  from  natural  soil  formation  or  through  technological  enhance¬ 
ment  of  the  natural  radiation  (23tl.l)>  Technologioal  enhancement  is 
responsible  for  elevated  indoor  radon  concentrations  in  homes  built 
o  or  reclaimed  phosphate  mina  landa  in  Florida  and  in  Colorado  where 
uranium  mine  tailings  were  used  as  fillers  ( 9 1 28  > .  Even  though  the 
general  populations  axpoaura  to  radon  progeny  has  been  increased  by 
mismanagement  of  mining  wastes,  the  overwhelming  majority  of  collective 
dose  equivalent  to  tha  us.  public  ia  due  to  the  natural  source  occur¬ 
rences  of  radon  and  its  progeny. 
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Typical  Ra-222  concentrations  in  U.S.  houses,  averaged  over  a 


year,  are  in  the  range  0.2  -  4.0  pCi/iiter  ( 24 1 27G ) .  This  value  has 
been  estimated  to  be  responsible  for  lung  dose-equivalents  that  are 
approximately  5  percent  of  the  oooupational  limit  ( 23 » i . 2 ) .  The  largest 
contribution  to  the  collective  lung  dose  equivalent  is  from  all  natural 
background  radiation. 

Prior  to  the  last  ten  years,  before  energy  conservation  measures 
were  instituted  in  home  oonstruotion,  typical  exposures  to  radon  progeny 
and  other  indoor  pollutants  were  low.  This  was  due  to  high  infiltration 
rates  usually  on  the  order  of  1  to  2.S  air  changes  per  hour.  8ince  that 
time  more  people  have  become  extremely  energy  eonseious  because  of 
increases  in  energy  oosts  and  fadarally-sponsorad  energy-conservation 
programs.  The  move  toward  more  energy-eff ioient  homes  has  led  to 
dramatio  reductions  in  residential  air  Infiltration  rates.  The  Depart¬ 
ment  of  Energy  estimates  that  over  50  percent  of  U.S.  home  owners  have 
taken  action  to  tighten  their  houses,  and  20-30  percent  have  an  air 
ventilation  rate  of  less  than  one  air  change  per  hour  (23 ■ 1.2).  The 
continuation  of  these  trends  and  the  fact  that  most  people  will  spend 
80  percent  or  more  of  the  day  indoors  (  39 1  277 )  will  mean  an  increase  in 
the  radon  progeny  concentration  and  greater  radiation  exposure  to  the 
general  public. 

Origin  of  Radon 

Radon-222  gas  and  its  short  lived  daughters  ere  the  intermediate 
decay  products  of  the  naturally-occurring  deoay  series  of  Uranium-238. 
The  parent  radionuclide  U-238  which  was  formed  at  the  earth's  creation 


is  relatively  stable  with  a  half-life  of  4.5  x  109  years  (  3  6 1  4  7 )  . 

Because  of  its  primordial  origins  tl-238  is  found  in  the  earth's  erustal 
material  in  oonoentrations  that  can  vary  greatly  with  location  and 
geological  formation  in  the  region. 

The  primary  radiation  exposure  from  the  Uranium  decay  series 
begins  with  the  decay  of  Radium-226.  Radon-222,  which  is  the  first 
decay  daughter  of  Ra-226,  is  an  inert  noble  gas  that  is  capable  of 
exhalating  from  the  air  spaces  within  the  soil  or  other  material  con¬ 
taining  the  parent  radionuclide.  Radon-222  has  a  half-life  of  3.82 
days  (36i43)  and  decays  into  the  elements  polonium,  lead  and  bismuth. 
The  primary  decay  scheme  for  radon  is  shown  in  Figure  1.  The  daughters 
are  also  referred  to  as  Ra-A,  Ra-B  and  Ra-C  respectively. 

Radon  exhalation  from  ground  soils  is  affected  by  several  physical 
and  meteorological  factors,  besides  the  soil  concentration.  The  condi¬ 
tion  of  the  soil,  porosity,  the  amount  of  air  spaces  found  between  soil 
constituents,  is  a  major  factor  in  determining  the  rate  at  which  gas 
moves  through  the  macerial.  To  illustrate,  the  diffusion,  movement  of 
the  gas,  through  sand  would  be  slower  than  the  diffusion  found  in  soil 
containing  large  fragments  with  large  interstitial  spaces.  The  moisture 
content  and  snow  or  ice  covering  on  soil  will  also  decrease  the  emana¬ 
tion  rate  of  the  radon  gas  ( 2 3 1 1 . 5 ) .  The  flux  (pCi/m*  *  sec)  is  reduced 
by  slowing  down  the  linear  diffusion  rata  which  causes  more  of  the  gas 
to  decay  before  it  escapes  into  the  free  air  space  above  the  ground. 
Because  the  radon  progeny  are  metal  atoms  they,  the  progeny,  will  then 
remain  trapped  in  the  medium.  Drastic  changes  in  the  tli.x  may  come 
about  as  a  result  of  sudden  transient  changes  in  barometric  pressures 
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above  tha  soil-air  interfaca  ( 1 7 i 2 1 6 ) .  A  drop  in  baromatric  pressure 
reduces  the  pressure  on  gas  emanation  and  exhalation  rates  can  increase. 
Under  this  condition  radon  gas  concentrations  could  be  several  orders 
or  magnitude  above  that  predicted  for  simple  diffusion  (23il.5).  Zn 
contrast  as  the  barometrio  pressure  increases  a  corresponding  decrease 
in  radon  gas  flux  could  result. 

Other  factors  such  as  temperature  and  winds  between  the  ground  and 
open  air  can  influence  the  flux;  however,  these  factors  are  of  leaser 
importance  because  they  may  occur  over  a  short  period  of  time. 

Sourosa  of  Indoor  Concentrations 

Radon  gas  exhalation  into  a  reduced  dilution  volume  will  cause  an 
increase  in  steady  state  concentrations  of  the  gas.  This  is  usually 
the  situation  when  homes  are  constructed  on  top  of  soil  containing 
radium  precursors  and  ventilation  rate  of  the  house  is  reduced.  The 
house  foundation  may  create  a  barrier  to  radon  gas  diffusion  causing 
more  of  the  gas  to  decay  before  it  moves  out  of  the  floor  or  medium. 
Radon  has  been  reported  as  having  a  relaxation  length  of  approximately 
8  cm  (Si  269  )  i  because  of  this,  large  length  radon  gas  easily  enters 
into  the  home.  The  soil  beneath  homes  is  generally  accepted  as  the  pri¬ 
mary  source  of  radon  entry,  However,  several  other  access  pathways  have 
been  found  to  be  important  sources  of  radon. 

Construction  Materials 

Construction  materials  which  contain  high  levels  of  radium  or 
uranium  contents  have  been  studied  and  reported  to  contribute  to  the 
indoor  concentrations  of  radon  gas  (13jl8»30).  These  studies  havo 


found  that  materials  auch  as  concrete,  camant  blocks,  gypsum  board, 
and  masonry  contain  radioactiva  matarial  which  contribute  to  tha 
source  of  radon  gas  in  buildings.  Concrata  made  from  beach  sands  and 
mixtures  of  granite,  pumice,  and  alum  shale  have  among  the  highest 
reported  radon  exhalation  rates  (26i307).  Another  source  of  high 
radium  content  in  building  materials  is  the  use  of  residues  from  indus¬ 
trial  processes.  A  large  contributor  ia  the  use  of  Phoaphogypsum,  a 
calcium  sulfate  produced  by  treatment  of  phosphate  ores  with  sulfuric 
said.  This  material  is  generally  used  in  wallboard.  In  this  treatment 
Ra-226  usually  follows  the  calcium,  leading  te  tans  of  pCi/g  in  the 
gypsum  (26i308).  The  U.S.  is  the  largest  producer  of  the  phosphate 
rook  material,  but  use  of  the  phosphogypsum  ia  found  in  other  countries 
having  little  natural  gypsum.  The  principal  use  of  a  byproduct  matarial 
in  the  U.S.  is  conarete  blocks  using  phosphate  slag,  whioh  contains  muoh 
of  the  Ra-226  and  U-238  in  the  phosphate  ore.  These  blocks  have  been 
estimated  to  have  Ra-226  concentrations  of  about  20  pCi/g  (26i308). 

Radon  in  Water 

Radon  gas  is  soluble  in  water,  but  the  solubility  decreases  with 
an  increase  in  temperature.  Ground  waters  which  are  encased  in  radium 
bearing  rock  formations  can  have  concentrations  of  Rn-222  greater  than 
30,000  pCi/1  (23 i 311).  The  radon  gas  can  quickly  transfer  to  the  air 
when  heated  or  when  aerated  as  in  a  shower,  a  faucet  aerator,  or  appli¬ 
ance.  Several  homes  in  Maine,  out  of  2000  surveyed,  were  found  to  have 
indoor  concentrations  of  radon  gas  in  excess  of  10  pCi/1,  with  the 
primary  source  of  the  radon  being  the  potable  water  (  12 1  34 4  ) . 
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Natural  Gas 


As  with  water,  daposita  of  natural  gat  located  in  or  near  soil 
formations  which  have  high  content  of  radium  can  also  become  a  source 
of  radon  gas  in  homes.  Radon  diffuses  from  the  rook  into  the  well  and 
infiltrates  the  gas.  Release  into  the  home  ooours  when  the  natural 
gas  is  used  for  cooking.  The  concentration  of  radon  gas  in  gas  distri¬ 
bution  lines  in  the  U.S.  has  been  found  to  range  from  1  to  100  pCi/1 
( 26 1 3 12 ) .  Differences  in  storage  and  transmission  times  will  alter  the 
radon  oonoentration  delivered  to  the  home. 

Indoor  Concentrations  of  Radon  and  Rrogany 

Ones  radon  gas  has  entered  the  home  it  will  mix  with  the  room  air 
and  decay  according  to  the  scheme  shown  in  Figure  1.  The  daughter 
produots  consist  of  heavy  metals  atoms  of  polonium,  lead  and  bismuth. 

Each  is  radioactive  and  possesses  a  short  half-life  less  than  30  minutes. 
These  decay  produots  attach  themselves  to  dust  particles  and  other 
aerosols  in  the  room  air.  Depending  on  the  dust  concentration  in  the 
room,  some  of  the  daughter  products  will  remain  unattaohed,  positively 
charged  ions  and  as  such  may  attach  themselves  to  other  small  polar 
molecules  like  water  vapor. 

Effect  of  Ventilation.  Many  homes  are  not  airtight  and  therefore 
outside  air  will  infiltrate  through  doors,  cracks,  and  windows.  This 
rate  is  governed  by  several  factors,  i.e.,  wind  speeds,  temperature, 
and  pressure  differences  between  inside  and  outside  air.  The  rate  of 
infiltration  air  into  the  home  has  a  great  bearing  on  the  radon  gas 
concentration  in  the  home.  In  a  home  with  a  volume,  V,  room  air  that 
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seeps  outdoors  at  a  rate  F  (volume/time)  must  be  replaced  by  an  equal 
volume  of  outside  air.  The  airborne  aotlvity  removed  during  a  time 
unit  per  unit  volume  is  then  7/V  *  (C(in)  -  C(out)),  where  C(in)  and 
C(out)  are  the  indoor  and  outdoor  oonoentrations.  Therefore  F/V  is 
the  removal  rate  oonstanh  or  room  infiltration  rate.  If  the  outdoor 
oonoentration  is  mueh  smaller  than  the  indoor  oonoentration  then  the 
removal  rate  beoomes  F/V  *  C(in)  <23il.9>,  which  implies  that  as  the 
infiltration  rate  inoreases,  the  radon  oonoentration  will  be  reduced  by 
the  faotor  F/V. 

Working  Level 

In  assessing  the  possible  insult  to  the  lungs  due  to  inhalation  of 
radon  progeny,  a  unit  known  as  the  Working  Level  (WL)  was  defined.  This 
is  the  link  between  progeny  oonoentration  and  the  eventual  alpha  energy 
imparted  to  the  lungs  per  volume  of  air.  One  Working  Level  is  defined 
as  any  combination  of  Ra~A,  Ra-B,  Ra-C,  and  Ra-C'  atoma/liter  of  air 
that  upon  daaay  will  release  1.3  x  10s  MeV  of  alpha  energy.  The  Working 
Level  in  a  room  containing  radon  decay  products  is  determined  asi 

WL  -  1.3  X  10*»Cl  ♦  5.07  x  10"*Ca  +  3.73  x  10_,Cj  ( 23 1 1 . 13 ) 

where  Cl#  Ca,  and  C,  are  the  concentrations  of  each  daughter  in  Poi/1 
The  potential  alpha  energy  in  MeV  imparted  to  the  lung  can  be 
calculated  fromi 

PAE(MeV)  -  1.3  x  10s  X  WL  ( 23  1 1 . 13 ) 

A  100  microcurie  source  j.  radon  in  equilibrium  with  its  daughters 
will  give  one  working  level  of  energy. 
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Effects  of  Air  Traatmant  Devices 

Control  or  reduction  of  radon  and  radon  progany  in  raaidanoas  can 
ba  grouped  into  thraa  broad  categories.  Tha  firat  of  thaaa  ia  to  raatrict 
tha  movement  of  radon  gaa  into  tha  homa.  Saoond  ia  to  provida  dilution 
of  tha  radon  gaa  oonoantration  by  inoraaaad  air  ventilation.  Tha  third 
oatagory  ia  removal  of  progany  through  tha  uaa  of  air  traatmant  davioaa 
which  in  addition  ramova  othar  harmful  air  pollutanta.  Thara  art  aavaral 
typaa  of  davioaa  whiah  onn  ba  uaad  in  tha  oontrol  of  radon  progany  and 
indoor  pollutanta.  Thaaa  includa  tha  alaotroatatic  praoipitator,  high 
afficianoy  filtar,  and  oailing  fan.  Reduction  in  radon  working  lavala 
uaing  auoh  davioaa  haa  baan  invaatigatad  by  aavaral  invaatigatora 
(2;  13i 14 ) .  Anothar  mathod  uaaful  in  tha  ramoval  of  radon  progany  that 
ia  receiving  ranawad  intaraat  ia  tha  uaa  of  ion  ganaratora. 

Soopa 

Starting  with  a  ataady-atata  oonoantration  of  radon  progany,  tha 
affaotivanaaa  of  an  alaotroatatic  praoipitator  ia  avaluatad.  Tha  modi- 
fiad  Taivoglou  mathod  ia  uaad  for  analyaia.  Only  limitad  oontrola  on 
humidity  and  tamparatura  wara  available*  tharafora,  axtanaiva  investiga- 
tion  aa  to  thair  contributiona  to  tha  ohangoa  in  oonoantration  lavala 
wara  not  maaaurad. 


general  approach 

Tha  approach  to  thia  invaatigation  inoludad  1)  tha  praparation  of 
an  experimental  chamber  for  obtaining  a  ataady-atata  aourca  of  radon 
gaa  and  radon  progany,  2)  refinement  end  verification  of  a  mathod  for 

i 

|  aampling  and  analyaing  radon  progany,  3)  collecting  of  radon  progany 

! 


under  varying  conditions,  e.g.,  with  and  without  air  treatment  devices, 

4)  analysis  of  results  under  varying  conditions  to  determine  effective¬ 
ness  of  air  treatment  device. 

Sequence  of  the  Report 

Chapter  II  contains  the  basic  theory  and  math  analysis  used  in 
this  thesis.  Chapter  ill  explains  the  radon  gas  generator  and  teat 
chamber.  Chapter  IV  contains  the  details  of  the  measurement  methodology 
omployed  and  the  equipment  used  for  the  experiments. .  Chapter  v  oontains 
a  raview  of  the  data  and  results.  Chapter  VI  contains  the  conclusions 
drawn  from  the  results  obtained  in  the  experiment  phase.  In  addition 
Chapter  VI  contains  a  list  of  recommendations  for  improvements  and 


further  studies. 


II.  Theory  and  Analysis 


Introduction 

Tha  daoay  product*  of  radon  gas  ara  formad  aa  poaitivaly  ohargad 
atom*.  Thaaa  atoma  raadily  attach  thamaalvaa  to  aaroaola  and  duat 
i  partiolaa  suspended  in  ambient  houaa  air.  Whan  air  ia  drawn  through  an 

air  sampling  davioa  tha  attached  and  unattached  daughters  are  deposited 
onto  a  filter.  The  filters  are  then  analysed  by  using  a  technique 
known  aa  gross  alpha  counting. 

Tha  concentrations  of  tha  short-lived  progeny  of  radon  attaahed  to 
dust  partiales  in  air  ia  measured  by  taking  three  integral  counts  of 
alpha  activity  oollaatad  on  the  filter.  In  a  technique  first  developed 
by  E.C.  Tsivoglou  ( 3 S « 4 0 )  and  later  modified  by  J.W.  Thomas  (34),  the 
concentrations  of  Ra-A,  Ra-B,  and  Ra-C  in  air  are  calculated  from  three 
integral  counts  of  the  combined  alpha  particle  activity  of  Ra-A  and 
Ra-C.  In  this  thesis  post  sample  count  intervals  of  2-5,  6-20  and 
21-30  minutes  were  used.  This  method  is  used  because  it  has  an  inherent 
advantage  in  that  a  simple  counter  can  be  used  to  reoord  activity,  aa 
described  later  in  Chapter  3. 

In  using  the  modified  Tsivoglou  method  the  following  assumptions 
were  madei  1)  only  radon  progeny  are  responsible  for  alpha  activity 
detected  on  filter,  2)  counter  efficiency  is  the  same  for  all  radionu¬ 
clides,  3)  concentrations  of  radon  progeny  remains  the  same  during 


sample  collection. 


Math  Methodology 


The  development  of  the  enalyaie  ueed  for  celeuletlng  the  progeny 
oonoentretlone  ie  accomplished  in  essentially  two  steps.  The  first  is 
to  solve  three  interrelated  linear  differential  equations  which  relate 
the  number  of  atoms  for  each  decay  product  deposited  on  a  filter  to  the 
airborne  radionuolide  atom  concentration. 

The  eet  of  differential  equations  which  describe  the  rate  at  which 

r 

radon  progeny  atoms  are  deposited  on  the  filter  isi 
Change  in  number  of  atoms  -  Production  -  Decay 


dnA ( t ) 
dt 


V  *  xi.i 


Vi111 


X.  ni(t> 


(1) 


where  i 


n 

X 

q 


i 

i 
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1  for  Ra-A 

2  tor  Ra-B 

3  for  Ra-C 

number  of  the  ith  type  of  atom  on  the  filter 
decay  aonstant  of  the  ith  type  of  atom  (min*1) 
air  concentration  of  the  ith  type  atom  (atom/liter) 
air  sampling  rate  (liter/min) 


The  equations  are  then  of  the  formi 


+  ay  •  f(x> 


(2) 


with  a  general  solution! 
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where  is  a  constant  of  integration.  In  the  solution  the  constants 

of  integrations  are  evaluated  using  the  initial  condition  that  each 
T* 


13 


,.C 


(t)  la  equal  to  zero  at  the  start  of  the  sampling  interval,  i.e., 

n°  (t)  »  0  at  t  •  0  .  In  a  sampling  tima  of  t^  tha  numbar  of  radc , 

progany  dapoaiead  on  tha  filtar  arai 
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Equation*  (4),  (5),  (6)  can  ba  writtan  in  matrix  form  a* 


or 
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where  capital  letters  in  aquation  (8)  denote  the  matrices  in  equation  (7). 
The  elements  of  the  matrix  X  are  given  in  Appendix  A. 

The  seoond  step  in  solving  for  the  airborne  concentration  is  to 
relate  the  expected  Ra-A  and  Ra-C  alpha  disintegrations  over  the  three 
counting  intervals  to  the  number  of  atoms  of  eaoh  radionuclide  present 
on  the  filter  at  the  end  of  sample  collection.  The  linear  differential 
equations  whioh  describe  this  relationship  are  also  of  the  form  in 
equation  (2),  with  a  solution  of  the  form  denoted  in  equation  (3).  rh ... 
evaluating  the  oonstants  of  integration,  the  time  t  is  taken  as  0  at 

!  . 

the  end  of  sample  collection,  and  the  number  of  atoms  on  the  filter  at 
start  of  counting  time,  n^(t),  is  set  equal  to  the  number  of  atoms  on 
the  filter  rt  the  end  of  sample  collection,  i.e.,  n^ ( t B )  *  rw(t)  at 
t  -  0  .  At  a  time  t,  measured  from  start  of  counting,  the  number  of 
eaah  type  of  radon  progeny  atom  on  the  filter  ia  given  by  the  relation¬ 
ships! 
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Taking  equations  (D)  and  (11)  and  integrating  ovt,r  the  counting 
intervals  £61  each  of  the  three  Integral  count  times,  one  obtains  total 
count  o£  alpha  particle  activity  on  the  filter  at  the  end  of  sample 
colleotion.  The  relationship  is  given  by  the  equation! 


^i)  "  t  \  !tl  (t)  dt 

L 


where  e  ■  produat  of  the  detector  end  filter  paper  efficiencies 

The  gross  counts  observed  in  the  detector  over  the  time  intervals 
are  related  to  the  expected  number  of  disintegrations  from  both  alpha 
emitters.  Substituting  equations  (9)  and  (11)  into  (12),  and  performing 
the  indicated  integration  results  in  an  equation  of  the  formt 

Cjtt)  -  e  nf(tg)  ♦  IjtCtj)  nj(t#)  ♦  lj5(t3)  nj(t§)  (13) 

The  total  counts  of  the  alpha  particle  aativity  on  the  filter  is  then 
related  to  the  number  of  radon  progeny  atoms  deposited  on  the  filter  at 
the  end  of  sample  collection  given  in  matrix  formi 
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The  matrix  elements  for  L  are  given  in  Appendix  A.  Combining  equations 
(8)  and  (15),  one  obtains  the  matrix  relationship) 


N  -  6.LVKQ 


(  16) 


The  solution  of  Aquation  (16)  for  Q,  tha  oonoantration  matrix,  results 
in  tha  matrix  aquation i 
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This  yields  a  aat  of  aquationa  relating  tha  three  integral  counts  of 
activity  on  tha  filter  to  tha  air  concentration  of  each  daughter.  The 
aquationa  arei 
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where  m  are  elements  of  the  matrix  M  derived  by  multiplying  the 
inverse  of  matrix  K  with  the  inverse  of  matrix  L. 


Analysis  of  Error 

In  any  experiment  there  is  introduced  in  the  final  result  a  certain 
degree  of  error  or  uncertainty.  These  uncertainties  can  be  attributed 
to  either  systematic  or  random  errors  in  the  experiment.  In  measuring 
the  airborne  concentrations  of  radon  daughters,  errors  are  introduced 
by  the  statistical  fluctuations  in  the  decay  process,  uncertainty  in 
the  regulation  of  air  flow,  and  the  efficiency  of  the  collection/detection 
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systam.  The  extent  of  uncertainty  in  the  radionuclide  atom  concentration 
was  calculated  through  application  of  the  standard  error  propagation 
formula  <  20 1 131 ) .  This  technique  assumes  that  the  counts  of  the  radio¬ 
nuclides  in  each  counting  interval  are  independent  of  each  other.  The 
error  analysis  here  includes  the  contribution  from  the  gross  counts  in 
each  of  the  three  time  intervals,  assuming  a  Poisson  distributed  counting 
error.  The  uncertainty  in  the  air  flow  rate  and  deteotor  efficiency  is 
also  evaluated.  The  error  of  the  calculated  concentration  of  the  radon 
progeny  is  then  estimated  by  the  equation  ( 19 1 12  )  i 


where  «  uncertainty  in  calculated  concentration  with  i  -  1  for  Ra-A, 
2  for  Ra-B,  3  for  Ra-C 

S  -  statistical  uncertainty  of  one  standard  deviation  in  counts 

®  3 

for  each  radionuclide  with  j  -  i,  2,  3 
Sv  ■  systematic  error  in  air  sample  rate 

>  systematic  error  in  detection  efficiency  of  counter 

The  systematic  uncertainties  in  the  air  sample  rate  and  detector  effi¬ 
ciency  were  calculated  to  be  4.5  and  4.9  percent  respectively. 


Weighted  Least  Square  Method 

A  weighted  least  squares  regression  analysis  for  calculating  the 


radon  progeny  concentrations  ot  each  radionuclide  deposited  on  a  filter 


by  an  air  sampling  device  was  developed  jointly  by  Dr.  D.  Shankland  and 


III.  Radon  Gas  Test  Chamber 


Introduction 

The  experiments  uonducted  during  the  course  of  this  investigation 
were  conducted  in  the  control  room  for  the  Hot  Cell  of  building  470 
Area  B,  Wright  ■'Patterson  APB,  Ohio.  Previous  air  samples  collected  in 
the  Hot  Cell  showed  concentration  levels  of  radon  progeny  to  be  in  the 
ranges  of  0.18-0.25,  0.16-0.2,  0.10-0. IS  pCi  per  liter  for  Ra-A,  Ra-B, 
and  Ra-C  respectively  (21i8a).  These  levels  were  too  low  to  effectively 
measure  any  reduotions  or  ahanges  in  concentrations  as  a  result  of  using 
air  treatment  devices »  therefore,  it  was  necessary  to  increase  the  radon 
and  radon  daughter  concentrations  by  artificial  meana.  The  purpose  of 
this  chapter  is  to  describe  the  method  and  equipment  used  to  enhance 
the  radon  gas  concentration.  Specific  areas  to  be  disoussed  are.i 

1)  Physical  Description  of  Room 

2)  Production  and  Distribution  of  Radon  Gas 

3 )  Environmental  Control 

Physical  Description 

The  radon  gas  test  chamber  was  set  up  using  a  procedure  similar  to 
that  described  by  S.N.  Rudnik  (29i464).  The  room  selected  was  the  per¬ 
sonnel  changing  and  service  room  of  the  Hot  Cell  in  Bldg.  470.  The 
physical  dimensions  are  shown  in  Figure  2. 

The  volume  of  the  room  was  calculated  to  be  99  m5.  The  room  had 
two  shower  stalls  and  two  toilet  stools,  with  all  inlet  water  lines 
sealed.  This  room  was  not  typical  of  most  living  areas  found  in  U.S. 


Dr.  G.  John  (30).  The  analyaia  program  waa  developed  to  calculate  the 
airborne  concent ration*  of  each  of  tha  daughter  products  and  will  be 
oompared  with  the  results  calculated  using  the  modified  Tsivoglou  method. 
A  comparison  of  the  daughter  concentrations  is  included  in  Appendix  £. 
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homes.  It  was,  however,  determined  that  tighter  environmental  controls 
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could  be  maintained  as  compared  to  other  available  rooms  in  the  building. 
The  walla  were  painted  cinder  block  and  the  floor  was  a  solid  concrete 
slab,  Physical  modifications  of  the  room  were  made  to  accommodate  the 
equipment  required  to  enhance  the  radon  gas  concentration.  The  major 
modifications  included  the  drilling  of  holes  in  the  wall,  for  the  purpose 
of  routing  necessary  air  and  powar  cables  to  the  equipment,  and  installa¬ 
tion  of  an  exhaust  duct  to  control  the  air  infiltration  rate  into  the 
room. 

Radon  Gas  Production  and  Distribution 

Radon-222  gas  waa  generated  from  a  100  mieroourie  souroe  of  Ra-226. 
The  radium  solution  was  made  by  mixing  a  100  microgram  liquid  source 
with  a  1  molar  solution  of  nitric  acid.  The  1  molar  acid  solution  was 
made  by  mixing  (9  ml  of  1.42  molar  solution  of  nitria  sold  with  1  liter 
of  deionised  water,  Ambient  air,  supplied  from  e  60  lb  air  source  was 
routed  through  a  water  separator  and  regulator.  The  air  supply  branched 
off  into  two  routes,  The  first  carried  air  through  the  radon  gas  system 
and  the  other  fed  directly  into  the  gas  manifold  providing  make-up 
(dilution)  air.  The  air  supply  lines  for  both  radon  gas  generation  and 
dilution  air  are  shown  in  Figure  3. 

The  air  supply  flow  rate  was  measured  using  a  Roger  Oilmont  Flow¬ 
meter,  size  2,  serial  number  B707.  ■  Tho  flow  of  air  through  the  radon 
gas  was  maintained  at  0.2  1/min.  The  air  was  passed  through  3  humidifiers 
before  bubbling  through  the  radium  solution.  The  humidifiers  were  mads 
of  1000  ml  flasks  filled  with  approximately  900  ml  of  water.  The 
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Figure  4.  Radon  Qaa  Ganarator 

introduoa  pra-humidif iad  air  to  tha  ganarator  and  tha  othar  uaad  to  trana- 
port  tha  radon  gaa  to  tha  diatribution  ayatam.  In  tha  radon  ganarator 
tha  humidifiad  air  flowad  through  a  final  humidifiar,  through  a  flaak 
containing  tha  radium  aolution,  through  a  liquid  oarry-ovar  trap,  and 
two  glaaa-wool-f illad  flaaka.  Tha  radon  gaa  waa  than  transportad  to  tha 
manifold  diatribution  ayatam  on  tha  floor.  Prahumidif ication  of  tha  air 
waa  dona  to  pravant  tha  radium  aolution  from  drying  and  contaminating 
tha  gaa  ayatam  with  radium  raaidua. 

Tha  manifold  ayatam  ahown  in  Figura  5  waa  conatructad  from  3.71  cm 
OD  poly-vinyl  chlorida  piping.  Holaa  with  0.0  mm  diamatar  wars  drillad 
in  tha  pipa  at  60  cm  intarvala.  Tha  pipaa  wara  placad  on  top  of  0.5 
inch  woodan  apacara  with  tha  holaa  facing  down  toward  tha  floor.  Tha 
ralativa  aiza  of  tha  hola  to  tha  diamatar  of  tha  pipa  inaurad  a  mora  avan 
diatribution  of  gaa  into  tha  room. 
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Figure  5.  Gas  Manifold  Distribution  System 
and  Ventilation  Exhaust  Duct 

Roughly  5  liter/min  of  dilution  air  was  added  direotly  to  the  mani¬ 
fold.  The  flow  rate  was  monitored  by  use  of  a  Roger  Qilmont  Flowmeter, 
size  3,  serial  number  F-116.  A  u-tube  Mercury  filled  manometer  was 
placed  in  series  with  the  radon  ga  i  faed  line.  The  primary  purpose  of 
the  manometer  was  to  monitor  a-y  blockage  or  restriction  of  the  air  flow 
through  the  gas  generator. 


The  radon  gae  coi.  lantration  was  moi.itored  uaing  a  3-liter  Lucaa  cell 
(Si 8-10 )  modified  Cot  'nuoue  air  flow.  T-u\  coll  was  *et  up  in  the 
middle  of  the  room,  with  t inlet  approximately  ;  cetera  above  the  floor. 

Environmental  Control 

For  a  ateady  rata  ef  radon  gaa  emanation,  the  indoor  concentration 
ia  highly  dependent  on  the  air  ventilation  rate.  Air  infiltration,  the 
uncontrolled  leakage  of  air  through  oraoka  and  gapa  in  the  building,  ia 
a  dominant  meohaniam  for  ventilating  houaea.  In  conventional  homea 
infiltration  ratea  vary  from  0.S  to  1.5  air  changea  per  hour  (24i278). 

The  relation  for  radon  gaa  to  vantilation  ia  given  by  the  following 
equation i 

«  (*?)  *  (*?)  * 

where  C  ■  concentration  in  exchange  air  from  outaida 
o 

e  •  emanation  rate  from  walla  and  floora 
S  -  area  of  aurfaoea  which  emanate  radon 
V  ■  volume  of  room 
Xv  »  air  exchange  rate 

XD  -  decay  conatant 

Aa  can  be  aeen  from  the  equation,  the  oonoentration  of  radon  gaa  ia 
inverae) /  proportional  to  the  air  ventilation  rate.  Becauae  of  the 
critical  dependence  upon  the  air  infiltration,  all  doora,  cracka,  and 
holea  were  aealed.  Additionally,  all  exiating  air  conditioning  and 
heating  venta  were  aealed.  To  inaura  that  a  ateady  atate  oonoentration 


TABLE  I 


Air  Sxhauat  Ratal  aa  Function  of  Damper  Foiition 


Damper 

Foaition 

Exhauat  Rata  (ACH) 

Condition  1 

Condition  2 

Condition  3 

0.0 

0.06 

0.03 

0.07 

0.5 

0.43 

0.39 

0.36 

1.0 

0.78 

0.76 

0.76 

l.S 

1.13 

1.17 

1.07 

2.0 

1.47 

1.43 

1.38 

2.3 

1.69 

1.71 

1.63 

3.0 

1.88 

1.88 

1.79 

3.3 

2.04 

2.07 

2.02 

4.0 

2.19 

2.27 

2.30 

3.0 

2.48 

2.66 

2.63 

of  radon  gaa  ii  maintainad  at  all  timea,  it  waa  imparativa  that  air  infil¬ 
tration  rataa  ba  tightly  controlled.  To  oontrol  tha  air  exchange  rata  in 
tha  room,  an  axhauat  fan  waa  inatallad  approximataly  13  am  abova  and 
oantarad  ovar  tha  antry  door,  aa  oan  ba  aaan  in  Figura  3.  Tha  fan  draw 
air  from  tha  room  and  axhauatad  it  to  tha  outaida.  Sinoa  thara  waa  no 
air  blowing  into  tha  room,  maka  up  air  had  to  aaap  into  tha  room  through 
tha  walla  and  availabla  apacaa  around  tha  door,  Tha  air  axhauat  rata 
from  tha  room  through  tha  axhauat  duct  waa  maaaurad  uaing  an  Anemaotherm 
aarial  numbered  7602.  Tha  Anamaotharm  waa  providad  by  tha  Baaa  Haalth 
Phyaica  aactioh.  Tha  maaauramanta  wara  mada  in  tha  duot  betwaan  tha 
chambar  wall  and  tha  axhauat  fan.  Tha  air  axchanga  rata  could  ba  varied 
from  0.06  to  2.46  air  changaa  par  hour. 

Tha  comp lata  air  axchanga  rataa  maaaurad  ara  liatad  in  Table  I. 

Tha  maaauramanta  wara  mada  under  three  different  condition!.  Condition  1 


was  without  tha  blower  on  the  electrostatic  precipitator  in  operation. 
The  other  two  were  with  the  blower  operating  at  low  apeed  and  high 
•peed.  As  can  be  seen  from  the  figures  in  Table  I,  not  much  difference 
in  the  flow  rate  was  introduced  by  the  blower.  The  damper  position  is 
a  measure  of  how  far  the  sliding  door  waa  opened  ( in  inohes ) . 


IV.  Measurement  Methods  and  Equipment 


Introduction 

Thia  chapter  oontaina  a  detailed  daaoriptlon  of  tha  mathod  and 
equipment  uaad  during  tha  oouraa  of  tha  experimental  phaaa  of  this 
atudy.  Tha  araaa  dlacuaaad  arai 

1)  Radon  Qaa  Concentration 

2)  Radon  Progeny  Concentration 

3)  Radon  Progeny  Removal  Davioa  Invaatigatad 

Tha  information  provided  oontaina  tha  axparlmantal  procedure  uaad 
to  obtain  the  data  collected,  and  tha  analyaia  technique  uaad  to  obtain 
tha  reported  raaulta.  Thua  tha  intent  of  thia  chapter  ia  to  provide 
auffioiant  detail  of  tha  maaaurement  methoda  ao  that  any  atudant  may  be 
able  to  repeat  tha  experiment. 

Radon  Qaa  Concentration 

Tha  staady-atata  oonoantration  of  the  radon  gaa  must  bo  accurately 
and  precisely  known  at  all  timaa  ainca  it  is  tha  parent  alamant  in  tha 
decay  series.  Any  statistical  variation  in  tha  concentration  will  di¬ 
rectly  add  to  tha  error  in  the  radon  progeny  which  already  has  error 
associated  with  counting  atatiatica. 

Sampling  Mathod .  To  monitor  the  radon  gas  concentration  within  tha 
teat  chamber,  three  separate  but  similar  methods  were  usedi  1)  Contin¬ 
uous  monitor  using  3.77  liter  Lucaa  Call  ( B 1 10-18)22 ) ,  2)  Grab  air  aample 
using  3.77  liter  Lucaa  Cell,  3)  Air  sample  using  a  Radon  Collection  and 
Transfer  System  (RCTS)  with  small  Lucas  Cell. 
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Continuous  Monitor.  In  the  continuous  monitor  msthod  a  large 
Lucas  call  was  placed  in  a  light  tight  canister.  The  Lucas  cell  was 
placed  directly  onto  a  4  inch  photomultiplier  tube,  which  was  connected 
to  a  set  of  counting  and  recording  equipment.  A  continuous  £low  of  air 
through  the  Lucas  cell  was  accomplished  by  connecting  the  Lucas  cell  to 
a  vacuum  line,  which  was  vented  above  the  roof  line  of  the  building.  A 
1/2  inch  oponing  located  on  the  top  of  the  cell  provided  an  inlet  for 
the  air.  To  insure  only  air  containing  radon  gas  was  entering  the  cell, 
a  paper  filter  was  placed  over  the  inlet  to  trap  any  radon  progeny  and 
dust  aerosols.  The  air  was  drawn  through  the  cell  at  a  rate  of  3  liters/ 
min.  The  air  flow  was  measured  using  a  site  3  Roger  Qilmont  Flowmeter. 

The  activity  of  the  radon  was  then  established  by  the  following 
process.  As  the  radon  gas  decays  within  the  Lucas  cell  the  alpha  emis¬ 
sions  cause  light  scintillations  when  they  interacted  with  the  Zn8(Ag) 
coating.  The  light  from  these  scintillations  are  detected  by  e  photo¬ 
multiplier  tube  to  produce  a  pulse  that  is  registered  as  a  count  by  the 
counting  equipment.  The  aounta  were  reoorded  for  a  1000  second  interval 
and  the  total  for  the  interval  was  printed  out  automatically  on  a  tele¬ 
type.  The  number  of  counts  detected  was  then  used  to  calculate  the 
radon  gas  concentration  in  the  test  chamber  using  an  integrated  count 
method  developed  by  Jonasaen  (  17*350)  .  This  method  does  not  require 
that  the  daughter  and  parent  radionuclides  be  in  equilibrium  to  calculate 
the  concentration  of  the  parent  radionuclide,  with  the  fast  flowrate  as 
compared  to  the  half-life  of  che  daughter  products,  equilibrium  would  not 
be  achieved  in  the  1000  second  count  period.  The  efficiency  of  the 
Lucas  cell  is  0.28  as  was  determined  from  data  collected  by  Qrosner 
( 3 1 12-16) .  Calculation  of  the  efficiency  is  in  'ppendix  u. 


Grab  Sample  with  Lucas  Cell ■  A  second  method  of  sampling  and 
monitoring  the  radon  gas  oonoantrationa  in  the  teat  chamber  involved 
uaing  another  large  Lueaa  cell.  The  procedure*  used  with  the  grab  sample 
were  quite  simple.  First  the  Luoas  cell  is  evacuated  with  a  vacuum  pump 
to  as  low  a  pressure  as  possible,  usually  200  microns.  The  evacuated 
cell  is  then  taken  to  the  test  chamber,  a  filter  was  placed  over  the 
inlet  valve,  again  to  insure  only  radon  gas  enters  the  cell.  The  inlet 
is  placed  through  a  50  mm  sample  port  in  the  aooess  door  to  the  test 
chamber.  The  inlet  valve  is  then  opened  allowing  the  radon  gas  to  enter 
the  cell.  It  is  oriticel  that  the  time  of  sample  collection  be  recorded 
as  this  will  be  needed  to  oalaulate  the  gas  concentration.  The  cell  ia 
then  plaoed  in  a  light  tight  oanister,  atop  a  photomultiplier  tube,  set 
up  similar  to  the  continuous  monitor.  The  oounts  are  then  recorded  for 
1000-aeoond  intervals.  The  samples  were  usually  oounted  for  a  period  of 
3-4  hours.  Again,  using  the  integrated  count  method  the  gas  concentra¬ 
tion  oan  be  determined. 

Radon  Collection  and  Transfer  System  ( RCT S ) .  The  third  method 
of  determining  the  radon  gas  concentration  was  used  as  a  verification  to 
tne  two  large  Lucas  cell  valuos.  The  method  incorporated  the  us*  of  the 
RCTS.  The  basic  idea  of  the  system  is  to  first  collaot  a  sample  t J  air 
from  the  test  chamber  ir.  an  evacuated  tank.  The  tic  is  then  sent  through 
a  series  of  traps  to  remove  the  water  vapor,  carbon  dioxide,  and  finally 
through  a  charcoal  filter  cooled  with  dry  ice  and  alcohol  where  the  radon 
gas  is  collected.  Once  the  radon  gan  has  been  separated,  the  charcoal 
is  heated  and  helium  ia  used  to  transfer  the  radon  released  to  a  Lucas 


csj.1.  The  Lucas  cell  is  than  pl&coa  on  top  of  a  50  mm  photomultiplier 
tube.  The  oounta  ace  then  recorded  every  1000  eeconde.  Thee*  counte 
ere  then  ueed  to  aelaulete  the  initial  gas  concentration.  A  more  detailed 
*tep-by-step  process  for  the  RCTS  is  in  Appendix  C.  A  computer  program 
which  uses  the  integrated  count  method  to  determine  the  gas  concentra¬ 
tions  is  listed  in  Appendix  D. 

Environmental  Air  Samples .  During  the  course  of  this  otudy 
measurements  of  radon  gas  were  made  using  the  equipment  and  methods 
described  Above.  The  measurements  were  made  in  an  attempt  to  insure 
that  no  significant  amounts  of  radon  gas  were  introduced  into  the  atmo¬ 
sphere  or  the  inhabited  shop  areas  in  and  arouhd  building  470.  Measure¬ 
ments  before  and  after  release  of  radon  gas  were  taken  at  various 
locations. 

Radon  Progeny  Concentration 

Sampling  Method.  Radon  progeny  concentration  in  the  test  chamber 
. arc  made  from  samples  aollected  through  e  sample  port  in  the  access 
door  leading  into  the  chamber.  The  sample  port  wee  simply  a  50  mm  hole 
cut  Into  the  plexiglass  window  in  the  door.  The  hole  wee  sealed  with 
air  conditioning  duct  tape  to  prevent  air  from  entering  tnc  room.  Radon 
progany  aamplas  wera  deposited  onto  a  50  mm  diamatar  membrane  millipore 
filter  using  a  Gelman  air  jamplar.  Tha  air  sampler  had  an  avtraga  air 
flow  rata  of  22  t  1  litar/min.  The  millipore  filter  has  essentially  a 
100  percent  collection  efficiency  for  attached  and  unattached  daughters, 
a.id  is  believed  not  to  contribute  signif icently  to  burial  losses  of 
alpha  emissions,  A  five  minute  sample  time  was  used  for  all  samples 
collected  during  this  thesis. 


Counting  Instrumentation .  Alpha  emissions  of  Ra-A  and  Ra-C  wets 
detected  by  alpha  scintillation  methods.  A  four  inch  photomultiplier 
tube  was  inverted  and  placed  on  top  of  a  ZnS(Ag)  phosphor  disc.  This 
apparatus  was  then  placed  in  a  light  tight  housing  over  a  sliding  tray. 
The  exposed  filter  paper  was  placed  on  this  tray.  The  collecting  side 
of  the  filter  was  than  moved  into  position  in  direct  contact  with  the 
znS(Ag)  phosphor  coating.  As  the  daughters  decayed  the  alpha  emission 
caused  scintillations  which  were  detected  by  the  photomultiplier  tube 
and  processed  through  the  detection  and  counting  equipment.  The  photo¬ 
multiplier  tube  was  operated  at  9S0  volts. 

The  efficiency  of  the  system  was  determined  using  a  22  mm  diameter 
RaD+£  diso  standard,  and  was  found  to  be  0.422  ±  .021  CPM/DPM.  Complete 
calculation  of  the  efficiency  is  in  Appendix  B.  A  Multichannel  Analyser 
made  by  Nuclear  Data  Incorporated  was  operated  in  the  multichannel  soal- 
ing  mode  to  monitor  and  record  the  counts.  The  counts  were  recorded 
for  20  second  intervals.  This  setup  was  useful  in  that  it  provided  a 
visual  indication  of  the  counts  detected  and  also  allowed  the  collected 
data  to  be  stored  on  0  incu  floppy  computer  disas. 

Counting  Method  and  Analysis.  The  calculation  of  the  radon  progeny 
with  the  gross  alpha  counting  technique  was  based  upon  the  Tsivoglou 
method  described  in  Chapter  II.  The  counting  inter' als  chooen  in  this 
study  were  the  2-3  minute,  6-20  minute,  and  21-30  minute  intervals  rec¬ 
ommended  by  Thomaa.  The  calculations  for  the  daughter  activity  concen¬ 
tration  were  programmed  in  M0ASIC  for  an  Epson  QX-10  computer.  This 
program  wil1  calculate  the  concantretiona  for  any  sampla  and  counting 
intarvals.  Tne  program  will  also  calculata  tha  uncartainty  in  tha 


concentrations  talcing  into  account  the  arror  in  air  flow,  detector  effi¬ 
ciency  and  counting  atatiatica.  Additionally  tha  program  will  calculate 
tha  value  for  tha  Working  Laval  baaed  upon  tha  daughter  concentrations . 

A  listing  of  tha  program  ia  provided  in  Appendix  D. 

Radon  Progeny  Removal  Device.  The  air  treatment  device  inveati- 

.1  ’  " 

gated  in  this  thesis  was  an  electrostatic  precipitator  or  electronic  eir 
cleaner  ( EAC ) .  The  EAC  used  in  this  study  is  a  commercial  produot  devel¬ 
oped  by  Honeywall  Corporation.  Zt  is  intended  to  be  used  as  part  of  the 
central  heating  and  air  conditioning  system  in  the  homo.  The  EAC  Con¬ 
sists. of  a  wire  mesh  pre-filter,  used  to  trap  large  partiole  and  debris,  : 
and  a  series  of  electriaally  charged  plates  where  the  radon  progeny  are 
collected.  The  entire  element  oan  be  removed  from  the  system  for  easy 
cleaning.  Tor  the  purpose  of  this  thesis  a  special  setup  was  built  by  a 
local  oontraator.  Tha  EAC  was  placed  in  a  duct  that  was  oonnscr.sd  di¬ 
rectly  to  s  high  volume  sir  blower  capable  of  drawing  air  through  at  a 
rate  of  34000  1/min.  Tha  electronic  eir  cleaner  reduces  the  radon  progeny 
concentration  by  creating  an  intense  electric  field  between  the  parallel 
plates.  As  the  bl> war  draws  tha  air  through,  the  duet  particles  pass 
through  the  plates  where  they  are  collected  and  thus  removed  from  tha  air. 
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Radon  Oaa  Conoantrationa 

Tha  anvironmantal  maaauramanta  for  radon  gaa  war*  m*da  uaing  tha 
3.77  litar  Lucaa  oall.  Tha  looationa  and  timaa  of  tha  maaauramanta 
bafora  and  aftar  atart  of  gaa  ganaration  ara  liatad  in  Tablaa  ZZ  and 
ZZZ  raapaotivaly .  Tha  tablaa  alao  liat  tha  valuaa  of  radon  gaa  conaan- 
tration  oalculatad  for  aaoh  aampla, 


TABLE  ZZ 

Ambiant  Radon  Oaa  Conoantration 
Lavala  Bafora  Ralaaaa  of  Radon  Oaa 


Looation 

Data/Tima 

Radon  Conoantration  (pCi/1) 

Workshop 

29  Oot/1818 

0.022  i  171 

Out a  Ida 

30  Qot/0844 

0.037  t  221 

Control  Room 

30  Oct/1718 

0.078  t  171 

Workahop 

30  Oct/1300 

0.023  i  9.1% 

TABLE  ZZZ 

Radon  Qaa  Conoantrationa 
Aftar  Start  of  Radon  Qaa  Flow 


Looation 

Data/Tima 

Radon  Conoantration  (pCi/l) 

Workahop 

7  Nov/1404 

0.0.70  ±  121 

Outaidi 

8  Nov/0951 

0.022  t  501 

Baaamant 

8  Nov/ 1610 

0.042  t  121 

Control  Room 

11  Nov/1720 

0.215  ±  54 

TABLE  IV 


Radon  Oaa  Conoantration  in  Taat 
Chambar  Bafora  Saaling  Air  Infiltration 


Data/Tima 

Radon  Conoantration  (pCi/L) 

5  Nov/1410 

0.611  ±  3.11 

6  Nov/0900 

1.96  *  4.01 

6  Nov/1431 

0.886  1  2.61 

7  Nov/0902 

1.41  ±  4.01 

9  Nov/1614 

0.647  *  3.44 

Aa  can  ba  seen  by  quiok  inapaotion,  tha  rala  «•  ohangaa  in  tha 
valuaa  in  Tablaa  II  and  III  ara  on  tha  order  of  0.01  pCi/litar.  With 
tha  axception  of  tha  valua  liatad  for  tha  oontrol  room,  all  wara  within 
tha  uncertainty  of  tha  valuaj  tharafora,  no  aignifioant  inoraaaaa  wara 
notad.  Tha  oontrol  room  ia  adjaoant  to  tha  radon  taat  chambar  and  for  a 
ahort  tima  air  waa  blowing  from  tha  taat  ohambar  into  tha  control  room. 

Tha  aouroa  of  tha  unwantad  air  waa  found  to  ba  air  blowing  into  tha  oham¬ 
bar  from  tha  raar  aarvica  antranoa  to  tha  building.  Thia  aouroa  waa 
aaalad  and  tha  blowing  into  tha  control  room  waa  atoppad.  Thia  additional 
air  alao  cauaad  larga  fluctuationa  in  tha  radon  taat  chambar  gaa  concan- 
trationa  making  it  impoaaibla  to  achiava  a  ataady  atata.  Valuaa  for  air 
aamplaa  takan  of  tha  taat  ohambar  ara  liatad  in  Tabla  IV.  Aa  ia  claarly 
aaan  from  Tabla  IV  ataady-atata  concantrationa  wara  not  achiavabla  bacauaa 
of  tha  influx  of  uncontrolled  air. 

During  tha  pariod  7-12  Dacambar  air  aamplaa  wars  again  takan  out- 
aida  tha  ahop  area.  Tha  aamplaa  wara  takan  uaing  tha  3.77  litar  larga 
Lucaa  call.  Tha  calculatad  values  for  tha  radon  gas  concentration  from 
tha  data  collacted  ara  liatad  in  Tabla  V. 
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TABLE  V 

Radon  Ota  Conoant rat  ion 
From  Outaida  Shop  Araa 


Data/Tima 

Radon  Concentration  (pCi/1) 

7  Dec/1420 

0.283 

± 

6.21 

10  Dec/0859 

0.448 

t 

14.91 

10  Dao/1427 

0.143 

t 

10.41 

11  Dao/0913 

0.113 

i 

13.01 

11  Dac/1538 

0.0008 

t 

15.01 

12  Dec/0856 

0.143 

t 

11.01 

12  Dec/ 1322 

0.104 

i 

13.01 

Tha  notad  inoraaaa  in  tha  gas  oonoantration  maaaurad  hara,  aa  oom- 
parad  to  valuaa  obtainad  in  October,  ia  a  result  of  changes  in  atmoaphario 
conditions,  primary  ohangaa  wara  in  outdoor  tamparatura  and  winds.  Tha 
ovarali  raault  of  thaaa  maaauramanta  ia  that  no  high  lavala  or  potantially 
dangaroua  laval  of  radon  gaa  waa  introduoad  into  tha  aurrounding  environ- 
mant  aithar  in  or  outdoora.  Tha  maximum  parmiaaibla  oonoantration  of 
radon  in  air  ia  3  x  10"*  microcurie/ liter  of  air  (6i490),  which  ia  3 
pCi/1 .  Tha  highaat  valua  maaaurad  with  our  Luoaa  call  waa  0.443  pCi/1, 
wall  within  tha  maximum  parmiaaibla  oonoantration  in  tha  air. 

Aa  waa  atatad  aarliar,  tha  maaauramant  of  tha  daughtar  concantrationa 
would  ba  mada  whan  a  ataady-atata  oonoantration  of  tha  radon  gaa  waa 
achiavad  in  the  chamber.  To  monitor  tha  gaa  concentration  in  tha  chamber, 
aeveral  air  aamplaa  wara  collected  with  tha  RCTS.  Tha  results  of  the 
aample  run  are  contained  in  Table  VI.  Tha  average  ataady  atata  valua  of 
tha  gaa  concentration  during  tha  couraa  of  thia  study  waa  9.98  ±  0.421 

pCi/litar . 
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TABLE  VI 


Radon  Oaa  Conoantrationa 
in  Taat  Chambar  Using  tha  RCT8 


Data 

Radon  Conoantration  (pCi/L) 

*14  Nov 

6.63  i  .131 

16  Nov 

12.19  ±  .166 

**21  Nov 

9.39  t  .186 

**29  Nov 

8.10  i  .201 

6  Deo 

10.23  t  .216 

Average 

9.98  ±  .426 

*  -  Samp la  oollactad  prior  to  aaaling  raar 
door  in  taat  chamber 

**  -  Sampla  collaotad  S  houra  after  atart  of 
air  traatmant  davioa 


Radon  Progeny  Maaauramanta  Blaotroatatio  Praoipitator 

During  tha  pariod  19  Novambar  ••  9  Daoambar  1904,  nina  diffarant  runa 
with  tha  electrostatic  praoipitator  wara  mada  to  measure  tha  effectiveness 
of  tha  praoipitator  in  raduoing  tha  working  laval  or  potantial  alpha 
anargy,  aa  a  function  of  radon  daughtar  conaantration .  Tha  runa  fall  into 
threa  oatagoriaai  1)  two  runa  mada  with  juat  tha  blower  fan  in  operation, 
2)  four  runa  mada  with  tha  praoipitator  and  blowar  in  operation,  3)  three 
runa  wara  mada  with  half  of  tha  air  intaka  to  tha  praoipitator  blocked 
off  in  an  effort  to  raduoa  tha  air  flow.  Tha  air  flow  through  tha  pre~ 
oipitator  waa  meaaured  to  be  approximately  34000  liter/minute  of  air. 

Thia  equatea  to  roughly  20  complete  air  ohangaa  par  hour. 

Category  1.  Tha  two  runa  in  thia  category  wara  mada  on  19  and  21 
Novambar.  In  each  of  thaaa  runa  only  tha  blowar  fan  waa  uaad  to  draw 


air  through  the  filter*  and  precipitator.  Initial  measurements  of  the 
daughter  concentrations  and  working  level  were  made  before  the  start  of 
the  blower.  The  fan  was  then  turned  on  and  allowed  to  run  for  IS  min¬ 
utes  before  the  daughter  samples  were  started.  The  samples  were  then 
oolleoted  approximately  every  thirty  minutes  thereafter.  Pigures  6  and 
7  show  graphically  the  reduction  noted  in  the  working  level  and  daughter 
concentration  as  a  result  of  the  fan  blowing.  The  reductions  are  eauaed 
by  a  trapping  of  the  daughters  in  the  pre-filter  of  the  precipitator  and 
to  ineraased  plateout  of  the  unatteohed  daughter r  on  the  wells.  The 
measured  reduction  of  both  samples  was  averaged  and  found  to  be  73  percent 
Category  £.  In  this  category  four  run*  were  made  to  measure  the 
effective  reduction  of  the  daughters  end  working  level  aa  a  result  of 
eleotrostatio  precipitator  operation.  The  first  run  was  made  on  26  No¬ 
vember.  As  with  the  measurements  in  oategory  1,  aamplas  were  first  col¬ 
lected  with  no  air  treatment  devices  in  operation.  The  SAC  was  started 
and  15  minutes  ware  allowed  to  pass  before  the  first  sample  with  the  air 
treatment  device  was  taken.  Subsequent  samples  at  thirty  minute  intervals 
wore  then  taken.  After  sample  8  was  collected  the  SAC  end  blower  fan 
were  turned  off.  One  hour  later  daughter  samples  were  again  taken  to  mea¬ 
sure  the  ingrowth  of  the  daughter  produots.  The  reduction  end  ingrowth 
for  this  run  is  shown  in  figure  8.  The  overall  measured  reduction  in 
working  level  was  90  peroent.  It  is  interesting  to  note  that  within  one 
and  a  half  hours  after  the  ZAC  was  stopped  the  working  level  and  daughter 
concentrations  were  almost  at  the  beginning  point  as  measured  before  the 
EAC  was  started.  The  average  reduction  for  each  daughter  was  SB,  93  and 
96  percent  for  Ra-A,  Ra-B  and  Ra-C  respectively. 


Figura  6.  Run  1  Plot  of  Working  Laval  and  Daughter 
Concantrationa  With  Blowar  Fan  Only 
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Figure  7.  Run  2  Plot  of  Working  Level  and  Daughter 
Concentration*  With  Blower  Fan  Only 
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Figure  e.  Run  1  Plot  of  Working  L«v*l  and  Daughter 
Concentratione  With  EAC  in  Operation 


Run  2,  made  on  2''  November,  was  identical  to  the  samples  collected 
in  Run  1;  the  only  difference  of  note  is  that  only  ona  sample  was  col¬ 
lected  with  no  air  traatmant  device  in  operation.  The  remainder  of  the 
procedures  were  exact!"  the  same.  As  shown  in  Figure  9  an  overall  reduc¬ 
tion  of  90  percent  in  working  level  was  experienced  aa  a  result  of  opera¬ 
ting  the  EAC.  The  corresponding  reductions  in  the  daughter  concentrations 
of  66,  9S  and  96  percent  for  each  of  Ra-A,  Ra-B  and  Ra-C  were  n.itad. 

Within  thirty  minutes  of  stopping  the  EAC  the  daughters  and  working 
level  had  again  raised  to  a  level  near  the  concentration  levol  at  the 
start  of  this  run. 

Kuna  3  and  4  were  made  with  the  same  proceduras  as  used  in  the  two 
previous  runs.  The  only  exception  is  that  run  4  samples  were  taken 
after  the  EAC  had  been  operating  for  a  period  of  24  hours.  This  was 
done  in  an  attempt  to  deteat  any  noticeable  degradation  in  the  effec¬ 
tiveness  of  the  pr joipitator  due  to  saturation  of  filter  or  long  opera¬ 
tions.  As  can  he  seen  in  Figures  10  end  11  no  significant  changes  were 
noted  in  the  operation  after  the  24  hour  period.  The  overall  reduction 
in  working  level  was  found  to  be  90  percent,  which  is  the  same  value 
found  for  the  two  previous  runs.  The  reductio  it  in  daughters  were  *>2 , 

92  and  97  percent. 

Category  3.  Ae  can  be  seen  in  Figures  8-10,  within  the  first  15 
minutes  of  operation  tha  EAC  had  reducad  tha  working  lavels  and  daughtar 
concentrations  to  tha  lowest  level  obtainable.  In  an  attempt  to  cry  to 
fit  tha  reduction  curve  within  the  first  15  minutes  the  air  intaka  was 
reduced.  This  waa  dona  by  sealing  one  hAlf  of  the  intake  with  paper. 

The  first  run  in  this  category  waa  made  using  the  same  sample  and  start 
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Figure  11.  Run  4  Plot  of  Working  Level  and  Daughter 
Concentrations  After  EAC  on  24  Houra 
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tirr.es  as  with  the  previous  samples,  i.e.,  15  minutes  after  start  of  EAC 
the  first  sample  was  collected.  Figure  12  shows  the  working  level  still 
reduced  to  the  lowest  obtainable  level  within  the  first  15  minutes.  The 
overall  reduction  measured  was  86  percent  for  the  working  level  and  50, 
77  and  97  percent  for  each  of  the  daughters.  For  the  next  two  runs  the 
time  delay  between  the  start  of  EAC  and  collection  of  the  first  sample 
was  shortened  to  5  minutes  with  subsequent  samples  taken  every  thirty 
minutes  thereafter.  From  Figures  13  and  14  a  gradual  redaction  in  the 
working  level  would  tend  to  indicate  the  reduction  is  linear.  The  over* 
all  reduction  in  working  levels  for  these  last  runs  were  93  and  90 
peraant  respectively.  The  average  reductions  in  the  daughters  for  these 
runs  were  84,  97  and  97  percent  for  run  2  and  67,  95  and  99  percent 
for  run  3. 

Daughter  Concentrations .  In  sach  of  tha  figurss  6  -  14  a  plot  of 
tha  radon  daughtara  concantrntion  versus  ths  Eleotronio  Air  Clerner 
oparation  time  ia  ahown.  The  variations  in  tha  daughtsr  oonoantration 
lavels  ara  a  rasult  of  the  high  oegree  of  disequilibrium  batwaan  the 
radon  gaa  and  its  daughtara.  Tha  raduction  in  Ra-A,  tha  uppar  lino, 
shows  an  average  reduction  of  56  percent  with  the  short  half  life  of 
3.05  minutes  and  30  minuta  daisy  between  sample  collection.  The  concen¬ 
tration  of  Ra-A  was  nearly  in  aquilibrium  with  the  radon  gaa  at  all 
timan.  The  second  two  drughters  with  half-lives  of  26  and  19  minutes 
are  never  in  equilibrium  with  their  parent,  as  tha  30  minute  sample  did 
not  allow  sufficient  time  for  ingrowth  of  the  daughters.  As  a  result  of 
this,  the  daughter  concentration  of  Ra-B  and  Ra-C  were  constantly  fluc¬ 
tuating  as  can  be  seen  in  each  figure. 
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Figure  12.  Fun  1  Plot  of  Working  Level  and  Daughter 
Concentrationa  With  Half  Intake  Sealed 
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VI.  Conclusion  and  Recommendations 

Conclusion 

The  objective  of  this  thesis  was  to  datarmina  tha  af feotiveneas  of 
using  an  alaotrostatio  praoipitator  as  a  maans  for  raduoing  tha  airborna 
oonoant rat  ions  of  radon  progeny  in  homas. 

To  aohiava  this  goal  savsral  praliminary  staps  had  to  ba  acoom- 
plishad.  Tha  first  of  thasa  waa  to  introduoa  radon  gas  into  a  tast 
ohambar  and  aohiava  a  ataady-stata  oonoantration  of  radon  gaa  and  its 
progany.  Dus  to  unoontrollad  air  infiltration  and  rsmoval  of  progany 
by  othar  msthods  bayond  tha  aoopa  of  this  pro j act,  a  quasi-ateady  atata 
waa  aohiavad.  Tha  avaraga  ataady  atata  gaa  oonoantration  during  tha 
oouraa  pf  this  projaot  was  9.98  pCi  par  litar  of  air  in  tha  ohambar. 

Tha  naxt  stap  waa  to  davalop  tha  thaory  and  analysis  to  oaloulata 
tha  daughtar  oonoantration  and  asaociatad  working  lavsl.  This  waa  dona 
uaing  tha  modifiad  Tsivoglou  msthod.  Tha  arrora  in  tha  calculated  valuaa 
wars  dua  to  tha  uncertainty  in  tha  air  flow  and  affioianoy  whioh  ware 
4.3  and  4.9  percent  raapaotivaly .  Also  errors  ware  introduced  by  tha 
counting  statistics  dua  to  randomness  of  tha  decay  process .  A  Poisson 
distribution  was  assumed  for  tha  counting  statistics.  As  a  result  of 
extrema  degrees  of  disequilibrium  in  daughter  concentration  tha  error 
associated  with  the  counting  statistics  was  notably  high.  Tha  errors 
in  tha  daughter  concentration  levels  ranged  from  12.8  to  431  percent 
for  Ra-A,  5.3  to  235  percent  for  Ra-B  and  9.8  to  531  percent  for  Ra-C. 

The  largest  error  in  tha  daughter  concentrations  was  noted  when  the 


disequilibrium  value  was  high  after  the  precipitator  had  been  on  for 
several  hours.  The  average  error  in  the  working  level  was  20  percent. 

The  final  step  in  meeting  the  objeotive  of  this  project  was  to 
take  air  samples  with  the  electrostatic  precipitator  in  operation  to 
measure  the  reduotion  in  the  daughter  concentrations  and  working  level 
due  to  the  precipitator.  In  all  of  the  measurements  that  ware  made  the 
precipitator  reduoed  the  working  level  by  as  muoh  as  90  percent.  This 
reduotion  usually  ooourred  within  19  minutes  after  start  of  the  precip¬ 
itator. 

From  the  date  oolleoted  and  analysed  results,  the  eleotrostatio 
precipitator  is  an  affeotive  method  to  use  in  reducing  the  radon  daughter 
oonoentrstions  in  homes. 

Recommendations 

In  an  effort  to  improve  the  results  and  better  evaluate  the  effec¬ 
tiveness  of  using  the  eleotrostatio  precipitator  the  following  recommen¬ 
dations  are  auggestedi 

1.  Submit  work  order  to  fabrication  shop  to  hava  seals  made  for 
all  heating  and  air  conditioning  ducts,  and  seal  off  ell  cracks  and  holes 
in  the  walla  and  doors.  Then  using  a  tracer  gas  make  detail  measurements 
of  the  air  infiltration  rates  in  the  room. 

2.  To  effectively  study  the  impact  of  humidity  and  temperature  on 
the  concentrations  of  radon  progeny,  submit  work  order  to  install  con¬ 
trols,  which  can  be  operated  from  outside  the  test  chamber,  which  con¬ 
trol  appropriate  equipment  to  regulate  the  humidity  and  temperature  in 
the  test  chamber,  i.e.,  dehumidifier  and  heaters.  A  humidifier  is  avail¬ 
able  with  the  present  setup. 
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3.  Determine  the  amount  of  daughters  which  are  removed  from  the 
air  due  to  increased  air  oiroulation  within  the  test  chamber  while  the 
blower  fan  is  in  operation.  This  will  give  a  more  accurate  value  of 
the  reduction  in  daughter  concentration  aausad  by  the  precipitator. 

4.  Using  the  present  setup  end  procedures,  continue  work  with  the 
weighted  least  square  regresaion  method  of  analysis  to  determine  the 
sensitivity  or  error  in  the  method,  Then  compare  the  error  results 
against  the  modified  Tsivoglou  method. 

5.  Conduct  studies  of  the  modified  Tsivoglou  method  under  condi¬ 
tions  of  extreme  disequilibrium  to  obtain  optimised  values  for  sample 
oolleetion  and  counting  times.  Conditions  of  extreme  disequilibrium 
will  be  achieved  within  an  hour  after  starting  the  precipitator. 
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Appendix  A i  Elements  of  the  K  and  L  Matrices  Used  in 
Calculating  Radon  Progeny  Concentrations 

This  appendix  lists  the  specific  expressions  for  the  elements  of 
the  K  and  1*  matrices  used  in  solving  for  the  airborne  concentrations  of 
each  radon  daughter  product. 

The  elements  of  the  K  matrix  used  in  calculating  the  number  of 


atoms  present  on  the  filter  after  sample  collection  aret 


Tha  elements  of  tha  L  matrix  uaad  to  calculate  the  number  of  alpha 


disintegrations  in  a  time  interval  tm  -  tnz.  Where  tm  and  tn2  are  the 
start  and  atop  count  times  for  each  post  sample  count  interval.  The 
L  elements  ares 
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Appendix  Bi 


Calculation  of  Efficiencies  for  the  Large 
Lucas  Call  and  the  Alpha  Counting  System 


The  calculation  of  the  large  lucaa  efficiencies  waa  made  using  a 
known  concentration  of  radon  gas  obtained  from  Monsanto  Mound  Labora¬ 
tories.  The  gas  oonoentration  was  SI. 74  pCi/liter  (8).  The  sample  was 
allowed  to  obtain  equilibrium  before  counting  started.  Tha  Lucas  cslla 
had  a  volume  of  3.77  liters,  giving  a  total  activity  of  193  pCi.  At 
equilibrium,  for  every  decay  of  radon  gee,  three  alpha  particles  are 
emitted. 

Lucas  cell  1  had  s  delay  of  20176  seconds  from  time  of  collection 
to  start  of  aountingi  the  activity  in  tha  call  was  than  187  pCi.  This 
gives  a  value  of  3.61  disintegrations  par  aaoond.  Lucaa  cell  3  had  a 
delay  of  73740  ascends,  giving  an  activity  of  167  pCi,  for  3.01  disinte¬ 
grations  per  second.  From  these  values  tha  sxpsntsd  counts  for  s  given 
time  interval  oan  be  determined.  Table  VII  lists  tha  expected  end 
observed  counts  for  s  20000  second  count  interval  and  the  efficiencies 
of  each  Lucaa  call. 

The  efficiency  of  the  groan  alpha  counting  system  'or  counting  the 
filter  papers  waa  determined  using  a  22mm  disc  source  of  RaD+E  which 
had  an  activity  of  308  dps  on  1  January  1933.  This  value  Is  on  file 
with  the  Physics  Laboratory  and  waa  providsd  by  tha  National  Bureau  of 
Standards.  Tha  expected  count  is  determined  from  the  aquation 

A  -  A0*  EXP  {-X  t)  .  Taking  X  to  be  8.3  x  10-sd_1  and  t  -  11634  daya. 
the  expected  activity  is  189  disintegrations  par  second.  For  a  count 
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TABLE  VII 


interval  of  4000  saconds  tha  axpaotad  counts  ar«  755880.  Tho  obsarvad 
oounta  for  tha  aourca  in  4000  saconda  was  319337,  r*»;ulting  in  an  abso- 
luta  affioianoy  of  0.422. 

Tha  arror  in  tha  affioianoiaa  can  ba  datarminad  by  application  of 
tha  itandard  arror  propagation  formula.  Oivan  a  5  parcant  arror  in  tha 
souroa  aotivitiaa  and  assuming  Poisson  diotributad  counting  statistics, 
tha  arrors  can  ba  aaloulatsd  quita  aasily. 


Appendix  Ci  Procedures  for  Use  o_f  the  Redon  Collection 

end  Transfer  Syetem  ( RCT3 ) 

This  appendix  will  explain  the  step  by  step  procedure  to  be  followed 
for  processing  an  air  sample  through  the  RCTS.  A  diagram  of  the  RCTS  is 
shown  in  Figure  15. 

The  purpose  of  the  RCTS  is  to  take  a  sampl*  of  air  assumed  to  con¬ 
tain  radon  gas,  process  the  sample  through  a  series  of  traps  to  remove 
the  water  vapor  and  carbon  dioxide  and  then  trap  the  radon  gaa  in  a  char¬ 
coal  filter.  Ones  reparation  of  theae  elements  has  been  accomplished, 
the  radon  gas  is  then  transferred  .into  a  Luoaa  Call,  where  the  radon  gea 
and  its  progeny  dauay  by  alph.  emission.  The  disintegrations  are  then 
detected  by  a  photomultiplier  tube  and  counted.  The  initiel  concentra¬ 
tion  of  the  gas  car.  thou  be  calculated  using  an  integrated  count  method 
developed  by  Jonaesen  (17i35C). 

To  inaura  the  radon  gae  collected  in  the  Luces  Cell  is  not  contam¬ 
inated  with  gas  from  the  ambient  air  in  the  RCTS,  it  ia  evacuated.  This 
is  done  by  opening  valves  2,  3,  8  and  allowing  the  vacuum  pump  to  draw 
all  the  Mr  from  tha  system  through  valve  8.  The  pump  ahould  be  allowed 
to  run  ur.ti-  the  vacuum  thormocouple  gauge  reads  about  30-40  microns. 

When  evacuation  is  complatsd,  close  valves  2  and  3. 

To  evacuate  th*  air  sample  collection  tank,  connect  the  tank  to  the 
RCTS  at  valve  1.  Open  valve  1  and  the  air  inlet  valve  on  the  air  tank, 
thin  will  allow  the  air  to  be  drawn  through  valve  8.  When  the  tank  is 
evacuated  then  close  all  valves.  The  next  stop  that  should  be  performed 
ia  to  evacuate  the  Lucas  Call  to  be  used  In  sampling  the  radon  gac .  This 
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Figure  13.  Diagram  of  tha  Radon  Collaotion  and  Tranafar  Syatam 

Is  dona  by  opanlng  tha  atopoook  on  tha  call  and  opaning  valva  S  on  tha 
RCTS.  This  faollltatas  tha  evacuation  through  v:uvo  5. 

After  tha  system.  collaotion,  and  sampla  ve***ls  have  bsen  evacuated, 
cold  bath*  ara  prepared  for  tha  water  and  radon  traps.  Tha  baths  ara 
preparad  by  mixing  mathyl,  athyl,  or  isopropyl  alcohol  and  dry  lea  in 
dewara.  Tha  trap*  should  ba  brought  dov/n  to  a  temperature  of  approxi- 
mataly  -7P’  C.  Tha  purposs  of  thasa  bath*  1*  to  frsaza  tha  watar  vapor 
in  tha  sampla  in  tha  watar  vapor  trap  and  to  fraaza  tha  radon  ga*  in  tha 
gaa  trap.  Tha  bath*  ara  slid  over  the  two  trap*  on  tha  RCTS  appropri¬ 
ately  marked.  Tha  time  required  to  bring  traps  to  tha  needed  temperature 
is  lpproximataly  one  hour. 
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The  tank  is  then  taker  to  the  location  where  the  air  is  to  bo  drawn, 


the  inlet  valve  is  openod  and  the  tank  is  allowed  to  fill  with  the  sam¬ 
ple.  The  tank  will  then  be  taken  to  the  RCTS  and  connected  at  valve  1. 

At  this  time  insure  all  valves  on  the  RCTS  are  closed.  Leaving  the 
inlet  valve  on  the  tank  closed,  open  valve  1  and  &  of  the  RCTS.  This 
will  evacuate  the  air  which  collected  in  the  feed  line,  and  could  contam¬ 
inate  the  sample  with  radon  gas  in  an  unknown  volume,  when  evacuation 
is  completed,  usually  in  45-60  seconds,  clots  all  valv*«. 

The  next  phase  of  the  process  will  bo  to  draw  the  eir  sample  from 
the  tank  and  pull  it  through  the  water  vapor,  carbon  dioxide,  md  radon 
traps.  Accomplishment  of  this  phase  is  done  in  the  following  steps. 

Open  valves  1.  2.  and  3  on  RCTS  and  the  inlet  valia  on  the  tank.  This 

will  pull  the  air  from  the  tank  up  through  the  radon  trap.  At  this  time 

gauge  B  will  zero  out.  The  purpose  now  is  to  pull  the  air  through  the 
traps  to  separata  the  radon  gas  from  the  other  elements.  This  ia  done 

by  regu.oting  the  air  flow  with  valve  4  which  is  connected  to  the  vacuum 

pump.  It  is  imperative  that  the  flow  through  the  system  ia  not  so  fast 
that  the  radon  gas  does  not  have  time  to  freeze  and  be  caugnt  in  the 
radon  trap.  It  recommended  that  a  minimum  time  of  15  mit.<  tea  for  34 
liters  of  aii  be  used  as  a  benchmark.  Again,  this  flew  is  regulated  by 
valve  4.  Wn«m  the  sample  has  been  drawn  through  the  traps,  the  valves 
ara  then  closed  in  the  following  sequence!  tank  inlet  valve,  vt,l<'e  1, 
valve  2,  valve  3,  and  finally  valve  4.  Record  the  time  r  .at  this  phase 
ia  complete  as  it  ia  critical  to  the  calculation  o£  the  gas  concentration . 
The  cold  b  th  is  now  removed  from  the  radon  chaico*.!  trap,  an  oven  is 
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raised  over  the  trap  and  voltage  applied  to  the  heating  element.  The 
temperature  is  raised  to  approximately  400*  C  which  takas  approximately 
60  minutes.  When  the  trap  has  reached  the  desired  temperature,  then 
the  next  phase  of  transferring  the  gas  to  the  Luoas  cell  will  begin. 

This  is  done  by  first  opening  the  valve  on  the  helium  tank.  Open  valve 
6  on  the  RCTS  until  gauge  A  reads  approximately  4  on  the  pressure  dial. 
Close  valve  6,  open  valve  7  which  allows  the  helium  gas  to  flow  into 
the  radon  trap.  The  helium  gaa  should  be  held  in  the  radon  trap  for  5 
minutes  before  transferring  it  to  the  Lucas  cell.  A  small  transfer  pump, 
which  pumps  the  gas  from  the  trap  into  the  Luoas  cell  is  started  at  this 
time  by  plugging  its  power  cable  into  an  electrical  outlet.  The  flow  of 
the  gas  to  the  Luoas  cell  is  regulated  through  the  TV  valve  identified 
on  sketch  in  Figure  IS.  While  leaving  valve  7  open  and  the  pump  running, 
slowly  open  the  TV  valve  until  the  needle  on  gauge  C  begins  to  move, 
leave  the  valve  open  until  the  needle  movement  stabilises.  At  stabili¬ 
sation  close  the  TV  valve  and  valve  7.  Open  valve  6  to  allow  more  helium 
gas  into  the  system,  repeat  the  process  until  gauge  C  reads  0,  indicating 
one  atmosphere  pressure  in  the  Luoas  cell.  When  transfer  has  been  com¬ 
pleted  the  Lucas  cell,  TV  valve,  valve  7  and  helium  tank  are  all  closed. 
At  this  time  the  Lucas  cell  is  removed  from  the  RCTS  and  placed  in  a 
photomultiplier  tube  and  counting  of  the  sample  can  be  started.  The 
counts  can  be  taken  for  any  length  of  time  desired.  The  gas  concentra¬ 
tion  is  than  calculated  using  an  integrated  count  method  described  by 
Jonaeaen  (  17*  350) .  The  efficiency  of  the  RCTS  has  besn  determined  by 
previous  experimenter*  and  is  taken  as  84  percent  ( 8  > 1 1  j  3 7  ) . 
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The  valua  of  84  percent  for  afficiancy  of  the  RCTS  waa  arrived  at 
by  procesaing  a  known  concentration  of  radon  gaa  through  tha  system. 

Tha  radon  gaa  waa  collected  by  flowing  air  through  a  known  concentration 
of  radium  in  solution.  Then  processing  the  gas  through  the  RCTS  using 
the  procedure  outlined  in  this  aypendix,  a  value  for  radon  gas  concen¬ 
tration  was  calculated.  This  calculated  value  waa  then  compared  to  the 
theoretical  value,  as  determined  from  known  radon  gaa  emanation  rate 
and  radium  solution  concentration.  The  value  of  the  efficiency  includes 
contribution  of  the  transfer,  collecting  and  counting  efficiency  of  the 
system;  therefore,  it  is  the  total  efficiency  of  the  system. 

In  processing  a  sample  of  air  with  an  unknown  radon  gaa  concentra¬ 
tion,  the  value  of  efficiency  for  the  system  may  be  higher  than  84  per¬ 
cent,  this  difference  resulting  from  the  faot  that  the  process  does 
not  include  the  efficiency  of  radon  gas  emanation  from  the  radium  solu¬ 
tion. 

It  is  recommended  that  the  efficiency  of  the  system  be  verified. 
This  can  be  accomplished  by  taking  an  air  sample  of  known  radon  gaa  con¬ 
centration  and  processing  it  through  the  RCTS.  The  sample  may  be  pro¬ 
vided  by  the  Monsanto  Mound  Laboratory  in  Centerville,  Ohio. 


Appendix  Di  Computer  Programs  Uaed  in  Hath  Analysis 


This  appendix  lists  the  computer  programs  used  in  performing  the 
calculations  for  this  thesis.  The  first  listing  is  the  program  uaed  in 
calculating  the  radon  gas  concentrations  using  an  intergated  count  method 
(16).  The  second  listing  is  the  program  for  the  modified  Taivoglou 
method  used  in  determining  the  daughter  concentrations.  Both  programs 
were  written  in  MBASIC  for  uae  on  an  Epson  QX-10  computer  system.  The 
third  listing  is  a  weighted  least  square  regression  method  used  in  deter¬ 
mining  the  daughter  concentrations.  It  is  writteh  in  Fortran  S  for  use 
on  a  Zenith  Z-100  computer  system  with  Z  DOS  operating  system. 


20  REM  THE  PROGRAM  USES  THE  INTEGRATED  COUNT  METHOD  FOR  THE  RADON  CON 
CENTRATION . 

30  REM  FOR  RADON  DETERMINATIONS,  THE  PROGRAM  REQUIRES.  1)  RADON  FREE 
LUCAS  CELL  BACKOROUNDCOUNT  AND  COUNT  INTERVAL  LENGTH  (IN  MIN)0t 
40  REM  2 )  TIME  (MIN)  BETWEEN  SAMPLE  COLLECTION  AND  COUNTING.  3)  COUNTING 
INTERVAL  (MIN),  AND  4)  TOTAL  COUNTS  IN  INTERVAL 
50  PRINT  "---LUCAS  CELL  DATA  INPUT-  —  " 

60  INPUT  "Luca*  background  interval  (in  Min)  «"|T2 
70  INPUT  "Luca*  background  count  -*>B2 
80  INPUT  "Lucaa  Call  Efficiency  -"»E2 
90  INPUT  "Tima  from  aample  to  count  (MIN)  ■*(T3 
100  INPUT  "Counting  interval  length  (min)  -*fT4 
110  INPUT  "Total  count*  -">G4 

120  Fr  23904. 80625#*(EXP(  - . 0001259*T3 )-EXP( - .  0001239*(T3+T4  ))  )-4 . 30748*(EX 
P( - .  2272*'  )-EXf <-.2272*(T3+T4)))-163. 51951#* ( EXP( -. 02386*T3 )-EXP(  - . 02588 
*( T3+T4 )  ) ) 

130  Q«F+93 . 656*  { i*XP(- .  03318*T3(  -EXP(  - .  03518* (T3+T4  ) ) ) 

3 40  R-(G4-B2*T4/T2)/(R2*G*2.22) 

150  D-100*(1/(2.22*«2*R*T4) )  *.  5 
160  Rl-R/34.41 
170  LPRINT  CHR8 ( 1 3 ) j 
180  PRINT  "  " 

190  PRINT  "—LUCAS  CELL  DATA—" 

200  PRINT  *T1  -  ",T3»*  Minutaa" 

210  PRINT  "T2-T1  ■  " » T4 » "  Minute*" 

220  PRINT  "Counter  Efficiency  -  *»E2f"  CPM/DPM" 

230  PRINT  "Total  Count*  •  "|G4 

240  PRINT  "Background  count  rat*  -  *>B2/T2» "+/-"|SQR(B2)/T2| "CPM" 

250  PRINT  "  " 

260  LPRINT  " 

270  PRINT  " — -Radon  Cone  and  Coeff  of  Variation---* 

280  PRINT  "Radon  -  "jRl»"PCI/L  +/-  "fD»"  *" 

290  LPRINT"  * 

300  END 
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10  REM  Program  to  calculate  the  radon  progeny  concentration  uaing  the 
20  RBM  Thomas  method  (Modified  Teivoglou  Method) 

30  DIM  A$(3), B|(3), A(  3 >  ,C(  3) , D(  3) ,  E{  3 ) ,  F(  3 ) , 0(3), Z{3) 

40  DIM  K(3,3).M3,3),M(3,3),S<3.3),R(3,3)IW<3).Y<3) 

50  DIM  B(2),T(6),U(3),V(3),P(3),J(3) 

50  INPUT  "Sample  oollaetion  tima"iT 

70  INPUT  "Air  flow  rate  (LPM)  and  uncertainty  in  l"rV,x 

60  INPUT  "The  affective  efficiency  of  the  ayatem( detector  and  filter) 

and  uncertainty  in  l"rQ,f 

90  PRINT  "The  counts  detected  for  each  time  interval" 

100  INPUT  "Counts  in  Tl"i<3<  1)  tPRINT 
110  INPUT  "Counts  in  T2  H»G(2) iPRINT 
120  INPUT  "Counts  in  T3  "i 0(3) (PRINT 
130  INPUT  "Background  count  rate  (CPM)  "»B 

140  PRINT  "Enter  the  start  time(tl)  and  end  tima(t2)  of  eaoh  count  interv 
el" 

150  INPUT  "T1  and  T2  for  first  interval" >T( 1 ), T( 2 ) 

160  INPUT  "T1  and  T2  for  seaond  count  interval  "jT(3),T(4) 

170  INPUT  ”T1  and  T2  for  third  count  interval  *»T(5),T(6) 

160  A*(l)  -  "Ra-A"  iBf(l)  -"first" 

190  A9 (  2 ) ■  "Ra-«V  iB*<2)  -"Ssoond" 

200  A$(3)  -  "Ra-C"  i B$ ( 3 )  -"third" 

210  D(l)-  .2272  t  D( 2  > -  .02586  i  D(3)»  .03518 
220  P(l)  -  D( 2 ) -D( 1 ) 

230  F( 2 )-  D ( 3 ) —  D ( 1 ) 

240  F(3)  ■  D(3)-D(2) 

250  FOR  I  -  1  TO  3 

260  C(I)  -  (Q(I)  -B*(T( 2*1) -T (2*1-1) ) ) 

270  Z(I)  -SQR( 0( I )  +  B* ( T( 2*1 )-T( 2*1-1 ) )  )/C(I) 

280  NEXT  I 

290  K(l,l)  -  ( 1-EXP ( — D ( 1)*T) )/D(l) 

300  K( 1 , 2  )  -  0 
310  K(l,3)  ■  0 

320  K( 2 , 1 )  •  1/D(2)-EXP(-D(1)*T)/(D(2)-D(1) ) +D( 1 )*EXP( -D( 2 ) *T } / ( D( 2 ) * ( D( 2 
)  -D(  1 )  )  ) 

330  K(  2 , 2  )  -  ( 1-EXP ( -D( 2 )*T) )/D(2) 

340  K( 2 , 3 )  -  0 

350  K( 3 , 1 )  -  (1-EXP(-D(3)*T) ) /D( 3 )-D( 2 ) *( EXP( -D( 1 )*T) -EXP( -D( 3 ) *T ) )/(F( 1) 
*F( 2 ) ) +D( 1 )*( EXP( -D( 2 ) *T )-EXP( -D(3)*T))/(F(1)*F(3)) 

360  K ( 3,2)  -  ( 1-EXP( -D( 3 ) *T ) )/D(3)  +  (  EXP)  -D( 3 ) *T )-EXP( -D(2)*T))/F(3) 

370  K ( 3,3)  -  ( 1-EXP ( -D( 3 )*T) ) / D ( 3 ) 

380  REM 

390  P.EM  Calculate  the  inverse  of  matrix  K 
400  GOSUD  1060 
410  FOR  I  -  1  TO  3 

420  L(I,1)  -  (1+D(2)*D(3)/(F(1)-F(2)))*(EXP(-D(1)*T (2*1-1) )-EXP( -D( 1 ) *T( 2 
•I) ) )+D( l)*D(3)/( F( 1)*F( 3) )*(EXP(-D(2)*T(2*I) )-EXP(-D( 2)*T(2*I-1) ) )+D( 1 )* 
D(2)/(F(2)*F{3) )*( EXP( -D(3)*T (2*1-1) ! -EXP ( -D( 3 ) *T( 2*1 ) ) ) 

4  30  L(I,2)  — D( 3 )/F< 3)*( EXP( -D{ 2)*T( 2*1-1) ) -EXP( —  D ( 2 ) *T( 2*1 ) ) )+D(2)/F( 3 ) * { 
EXP( -D( 3 ) *T( 2*1 ) ) -EXP( -D( 3 )*T( 2*1-1 ) )) 

440  MI.3)  -  (EXP(  -  D(3)*T  (2*1-1)  )-EXP(-D(  3  )*T(  2*X  )  )  ) 
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450  NEXT  I 

460  REM  Calculate  tha  inverse  of  matrix  L 
470  OOSUE  1310 
480  RKM 

490  REN  Multiply  matrix  a  with  matrix  r 
500  REM 

510  FOR  I  -  1  TO  3 

520  FOR  J  ■  1  TO  3 

530  M(I,J)  -  0 

520  FOR  K  -  1  TO  3 

550  M ( X , J )  -  M( I , J )  ♦  S( I , K) *  R( K, J ) 

560  NEXT  K 
570  NEXT  J 
580  NEXT  I 
590  REM 

600  REM  return  matrix  m  ■  s*r 
610  REM 

620  REM  multiply  matrix  m  with  a 
630  REM 

640  FOR  I  -  1  TO  3 

650  V ( I )  >  0 

660  FOR  K  -  1  TO  3 

670  Y(I)  ■  Y( I )  +  M(I,K)  *  C(K) 

680  NEXT  K 

690  NEXT  Z  ' 

700  REM 
710  REM 

720  REM  multiply  matrix  y  with  1/gv 
730  REM 

740  FOR  I  •  1  TO  3 

750  W( I )■  (1/(0*V) )*Y(I) 

760  NEXT  I 
770  REM 

780  FOR  I  -  1  TO  3 

790  J(I)  -  ((Z(I)*C(I))/<0*V)r2 

800  NEXT  I 

810  FOR  I  -  1  TO  3 

820  U(I)  -  SQR(M(I,ir2(*J(l)  +  (M(I,2r2*J(2)  +  (M(I,3r2)*J(3)) 
830  V(I)  -(U(I)*100)/W(I) 

840  P(I)  -  SQR(  V(  I  )A2+Y*2  )  ♦X,'2  )  ) 

850  A ( I )  -  (D(I)*W(I) )/2.22 
860  NEXT  I 
870  E< 1)-  11 

880  E<  2  )  -  (  D(  2  )  *W(  2  )  )  /  (  D(  1 )  *W  (  1 )  ) 

890  E(3)  -  ( D( 3 ) *W ( 3 ) ) / ( D( 1 ) *W ( 1 )  ) 

900  REM 

920  REM  print  solutions 
920  REM 

930  FOR  1  -  1  TO  3 

940  PRINT  “*****************  RADON  PROGENY  CONCENTRATIONS  *** 

**» 

950  PRINT  "Counts  for  ";B$(Ib"  interval  are  ">  G(  I ) 


960  PRINT  "Concentration  is  ";A$(I)>  "  ia  "  >  A(  I )  i  "pCi/llter  +/-  "jP(I) 
970  PRINT  "Equilibrium  concentration  between  Ra-A  and  “»A$(I )»"  ia"E(I) 

980  PRINT  "«*»*»»**«**«**»«*#***»*«*******»#***♦**«**##♦**•**»*»***** 
♦***####** 

990  PRINT) PRINT 
1000  NEXT  I 

1010  PRINT  "a*********************  WORKING  LEVELS  a*********************' 
***** 

1020  WI  -  .001025«*A(1)+.0030624*A(2)+.0037247*A(3) 

1030  SI  -  SQR( (. 0010256 )A2*( P( 1  )*A{ 1) /ICO ) A2+( , 005624 ) A2*( P( 2 )#A(2 )/  100)A 
2+(.0037247)A2*(P(3)*A{3)/100)A2) 

1040  PRINT  "Working  level  «  "»WI»M  +/-  "j  (SI/WI )*100» *4" 

1090  END 

1060  FOR  I  -  1  TO  3 
1070  FOR  J  -  t  TO  3 
1080  3(1, J)  -  K( I , J ) 

1090  NEXT  J 

1100  NEXT  X 
1110  FOR  I  ■  1  TO  3 
1120  N  -  3(1,1) 

1130  FOR  J  -  1  TO  3 

1140  IF  J  •  X  THEN  1170 

1150  3(1, J)  -  8{ I, J )/N 

1160  GOTO  1180 

1170  0(1, J)  -  1/N 

1180  NEXT  J 

1190  FOR  J  -  1  TO  3 

1200  IF  J  m  I  THEN  1280 

1210  Q  -  -S(J,I) 

1220  FOR  K  *  1  TO  3 

1230  IF  K  ■  l  THEN  1260 

1240  S(J,K)  -8<I,K)*Q  +  8(0, K) 

1230  GOTO  1270 

1260  S(J.K)  -  S ( I , K )  *  Q 

1270  NEXT  K 

1280  NEXT  J 

1290  NEXT  I 

1300  RETURN 

1310  FOR  I  -  1  TO  3 

1320  FOR  J  -  1  TO  3 

1330  R(I,J)  -  L( I , J  ) 

1340  NEXT  J 
1350  NEXT  I 
1360  FOR  I  -  1  TO  3 
1370  N  •  R( I , I ) 

1380  FOR  J  ■  1  TO  3 
1390  IF  J  -  I  THEN,  1420 

1400  R(I,J)  -  R( I, J  )  /N 


1410  GOTO  1430 
1420  R{ I , J )  -  1/N 

1430  NEXT  J 

1440  FOR  J  •  1  TO  3 

1450  IF  J  ■  I  THEN  1530 
1460  Q  -  -R( J, I ) 

1470  FOR  K  -  1  TO  3 

1480  IF  X  -  Z  THEN  1510 

1490  R(J,K)  -  R(I,K)  *  Q 

1500  GOTO  1520 

1510  R ( >T , K )  -  R(I,K)  *  Q 

1520  NEXT  K 

1530  NEXT  J 

1540  NEXT  1 

1550  RETURN 


non  no  non 


PROGRAM  RADON  MAIN 

SOLVES  FOR  THE  INITIAL  DISTRIBUTION  OF  A  RADIOACTIVE  CHAIN 
GIVEN  SEQUENTIAL  VALUES  OF  COUNTS  EACH  FOR  THE  SAME  TIME 
INTERVAL. 

DIMENSION  A(10)  ,Y(10)  ,B(10)  ,W(10)  fU(10,10)  ,ZZ(10,10), 
DIMENSION  Q {10}  ,CC(10,10)  ,DA(10)  ,LC(10)  ,UU(10)  ,TQ ( 10) 
DIMENSION  YY(96) ,AL(10) ,P(10) 

CHARACTER* 6 4  FNAME 

WRITE  (*,' (A\) ’) •  INPUT  FILE  NAME? 

READ {*,’ (A) V) FNAME 
NA-10 

Nt  *  OF  NUCLIDES  IN  THE  CHAIN 
A:  DECAY  CONSTANTS  OF  EACH  DAUGHTER 
W>  *  OF  ALPHAS  PER  DECAYING  NUCLEUS 
N-3 

A(l)-0. 0037877 
A  ( 2) ■0 .00043106 
A(3)-0. 00058642 
W(l)  ■  1 
W(2) *  0 
W  ( 3 )  -  1 
ST-3 

CONSTRUCT  EIGENVECTORS  OF  THE  DECAY  MATRIX 
*********************.«  ************************************ 

CALL  EIG  (U,A,NA,N) 

DO201I-1 *N 

WRITE (* , 899)  (U(I, J) , J-1,N) 

********************************************************** 

CONSTRUCT  TRANSFORMED  WEIGHTS  AND  ZERO  SOME  ARRAYS 
********************************************************** 

D014I-1,N-1 
D014 J-I+l , N 
W(I)-W(I)+W(J) 

DOl  I-1,N 

Y  ( I ) -VPROD (U  ( 1 , 1 )  ,1,W,1,N,XX) 

WRITE  (*  f  898)  (Y(I),I-1,N) 

C-0.0 
D013I-1,N 
Q  ( I ) "0 . 0 
W ( I ) -0.0 
B(I) -0.0 
D013J-1,N 
ZZ (I , J) -0,0 


READ  DATA 


WRITE ("*991) 
READ {*  *  992) F 
WRITE (*  *  993) 
READ (* , 994) EFF 


nnon 


WRITE  (*,996) F 
(C-v  WRITE  (*,997)  EFF 

''■<  Ml  #  OF  INTERVALS  AND.  COUNT  DATA 

WRITE (* , 902)  , 

READ  (*  ,  904) M 
WRITE  (*,906)  M 

OPEN ( 7, FILE-FNAME, ACCESS- ' SEQUENTIAL' , STATUS" ’ OLD ’ ) 
READ(7, 903 , END-40) ( YY ( I) ,1-1, M) 

40  CONTINUE 

CLOSE ( 7 , STATUS- ’  KEEP ' ) 

C  TNI  TIME  OF  START  OF  COUNTING  FROM  SAMPLE  COLLECTION 

WRITE  (*,935) 

WRITE  (*  ,  945 ) 

WRITE  (*,955) 

READ  (*,900)  TN,DT 
WRITE  (*,965)  TN ,DT 
WRITE  (*,975) 

VALUES  OF  "TM  AND  "D"  NOW  ENTERING 
Tl  TIME  AT  END  OF  COLLECTION  OF  "D"  COUNTS 
DS  COUNTS  ACCUMULATED  IN  EACH  INTERVAL 

T-TN 

WRITE (*,976) 

D03  K-1,M 
T-T+DT 

WRITE  (*,985) K 
D-YY(K) 

WRITE  (*  ,  995)  T,0 

C  CONSTRUCT  "DECAYED  WEIGHT  VECTOR" 

D04  I-l ,N 

4  Z ( I ) » (EXP ( -TN*A ( I) ) -EXP (-T* A( I) )  ) * Y ( I ) 

TN-T 

C  ACCUMULATE  RIGHT-HAND-SIDS  AND  COEFFICIENT  MATRIX 

C-C+D*D 
WRITE (* , 921) C 
D05  I -I , N 
B ( I ) -B ( I ) +  D*Z ( I) 

W ( I ) -B  ( I ) 

D05  J-l ,  N 

5  ZZ(I,J)«ZZ(I,J)+Z(I)*Z(J) 

3  CONTI  t/Ji 

WRITE  (*,924) 

DO  3. 01  -1 ,  N 

10  WRITE (*,925) (ZZ(I,J) ,J-1,N) ,B(I) 

C  SOLVE  FOR  INTERMEDIATE  VECTOR 

DO20I -1 ,  N 
DA (I) -0.0 
DO20J-1 ,  N 

20  CC(I,J)—  U(J,I) 

CALL  QU AD ( B , ZZ , CC , DA , Q , UU , N , N , N A , LC , TQ ) 


o  o  o 


C  COMPUTE  SENSITIVITY 

Q  Aft**************************************************** 

C-0.5* (C-VPROD  (W,1,Q,1,N,XX) ) 

WRITS (*,922)C 

C  FINISH  INVERSE  OF  ZZ 

ZZ(1,1)*1./ZZ(1,2) 

DO  11  J  -  2,N 
ZZ(J,J)*1./ZZ(J,J+1) 

DO  II  I-  J-1,1,-1 

11  ZZ  (J, I) "-VPROD (ZZ (I , J+l) ,-NA, ZZ ( J, J) ,-NA, J- I , XX) *ZZ  (I  ,  I) 
NP-N+1 

DO  12  J*1,N 
DO  12  I-J,N 

ZZ(I, J)«C*VPROD(ZZ(I,J) ,1,ZZ(I,I) ,1,NP-I,XX) 

12  ZZ(J,I)-ZZ;i,J) 

*******************************w********************** 

UNFOLD  INTERMEDIATE  VECTOR  TO  FINAL  SOLUTION 

ft***************************************************** 

D06  1*1, N 
X-0.0 
D07  J»l, N 

7  X-X+U(X, J)*Q;J) 

Z(I)*X 

AL(T)-Z(I)/(F*EFF*ST) 

P(I)-(AL(I)*A(I) )/0.037 
>'  6  CONTINUE 

C  PRINT  RESULTS 
WRITE (* , 920) 

WRITE (* , 930) N, (Z(I) ,I*1,N) 

WRITE (* , 960) 

WRITE (*, 970) N, ( AL ( I ) , I *1 , N) 

WRITE-C1 , 980) 

WRITE(*, 970)N, (P (I) ,1-1, N) 

WRITE  (*, 940) 

DO  2  1*1, N 

2  WRITE(*, 950)  (ZZ  < I , J) ,J-1,I) 

STOP 

99  FORMAT (I5/(2F15. 5) ) 

98  FORMAT ( 1 5 ,F10 . 0) 

97  FORMAT (2F10.0) 

899  FORMAT ( 1H  ,'U  -  ',  G14.7) 

898  FORMAT ( 1H  ,'Y  -  ',  G14.7) 

900  FORMAT ( 2F5 . 0 ) 

903  FORMAT (10F6.0) 

'991  FORMAT (1H  ,' ENTER  FLOW  RATE  IN  LITER/MI N  IN  F5.0  FORMAT') 

992  FORMAT (F5 , 0) 

993  FORMAT  (1H  ,' ENTER  THE  PRODUCT  OF  DETECTOR  AND  FILTER  EFf) 

994  FORMAT (F4. 3) 

996  FORMAT ( 1H  ,'FLOW  *',F5.0) 


W-JN.  W  -*  ..Hi  r*-*-r'*  - 


FORMAT (1H  , 

FORMAT (1H, ' 

FORMAT (14) 

FORMAT (F7.0 
FORMAT (1H  , 

FORMAT (1H  , 

FORMAT (1H  TN-  TIMS  FROM  END  OF  SAMPLING  TO  START  OF  COUNTING') 
FORMAT (1H  DT-  COUNTING  TIME  PER  MEASUREMENT  IN  SECONDS') 

FORMAT (XH  ,  '  TN,DT-  ',2F9.0) 

FORMAT (1H  VALUES  OF  "T"  AND  "D"  NOW  ENTERING’) 

FORMAT ( 1H  ,  '  ENTER  JUST  D,  COUNTS  AT  TIME  T,  IN  FREE  FORMAT') 

FORMAT ( 1H  ,  '  TIME  AND  COUNT  SET,# - ',14) 

FORMAT (1H  ,  '  T,D  ■  ' ,2F9.0) 

FORMAT (1H  , 'ZZ  -  MATHIX,  B  -  VECTOR'  ) 

FORMAT (//2G14. 7) 

FORMAT (1H  ,  'C  ■  '  ,F14.0) 

FORMAT (F14. 6) 

FORMAT (1H  ,' NUMBER  OF  ATOMS  OF  EACH  RADIONUCLIDE') 

FORMAT (//I4/(6F12.0) ) 

FORMAT (1H  , 'CONCENTRATION  IN  ATOM3/LITER' ) 

FORMAT (//I 4/ (6F12 . 4 ) ) 

FORMAT ( 1H  , 'CONCENTRATION  IN  PCI/LITER’) 

FORMAT (1H  ,'  MEAN  SQ.  ERROR  MATRIX') 

FORMAT (2X,1P6E13. 4) 

END 

SUBROUTINE  EIG (U , A, NA , N) 

CONSTUCTS  THE  EIGENVECTORS  OF  THE  COUPLED-DECAY  MATRIX 
DIMENSION  U (NA,NA) , A (NA) 

ZEROING  THE  ARRAY 
DOl  I-1,N 
DOl  J*1 , N 
U(I,J)-0.0 

STEPPING  THROUGH  THE  COLUMNS  (EIGENVECTORS) 

D02  J-1,N 

FILLING  THE  BOTTOM  PART  OF  THE  VECTOR 

D03  I-J,N 

X-1.3 

IM-I-1 

IF  ( IM-J )  6,7,7 
DO 4  K*J , IM 
X-X*A(K) 

IP-I+1 

IF  (N-IP) 3,8,8 
DOS  K-IP,N 
X-X* (A(K)-A(J)  ) 

U  ( 1 ,7)  "X 
CONTINUE 

NORMALIZE  VECTORS 
D09J"1 , N 


, 'EFFECTIVE  EFFICIENCY  -',F4.j) 

'  ENTER"M" ,  NO  OF  DATA  IN  "14"  FORMAT'/) 


'M“ ' ,15) 
'ENTER 


"2F5.0"  FORMAT') 


'v/77 


DO10I-1,N 

X-X+U(I,J)**2 

X-1./SQRT(X) 

D09I-1,N 

o(i,j)-x«o(i,cn 

RETURN 

END 

SUBROUTINE  QUAD (A,B,C,D,X,U,NX,NC, NB,LC,TN) 

DIMENSION  A(l) ,B(NB,1) ,C(NB,1) ,D(1) ,X<1) ,0(1) ,LC(1) ,TN(1) 
NP-NC+1 

REDUCE  THE  PROBLEM  TO  "NORMAL"  FORM 
CALL  MCDCMP(B,NB,NX,X) 

D01  1*1 » NX 
C (I |NP) *A(I) 

D02J*1 , I 
B ( Jf 1+1) *B (I f J) 

DC3J-1,NP 

C(1,J)*(C(I,J) -VPROD (B  ( 1 , 1 )  , NB ,C ( 1 , J) , 1 , 1-1 ,XX) ) /B ( I , I ) 

CONTINUE 

D07J-1,NC 

LC ( J) *J 

U  ( J) -0.0 

XX-l./SQRT (VPROD (C(1,J)  ,  1  ,C  ( 1  ,  J)  ,  1 ,  NX  ,  XX) ) 

D09I *1 , NX 
C(I,J)-XX*C(I,J) 

D(J)-XX*r  (J) -VPROD (C(l, J) ,1,C(1,NP) ,1,NX,XY) 

TN ( J ) -XX 
CONTINUE 
NP-1 
NI  "0 

MAIN  ITERATION 
D04 1 *1 , NX 

X (I ) ■  VPROD(C (I ,1) ,NX,U,1,NI ,XY) 

PREPARE  FOR  LM  COMPUTATION 
00281*1, NI 

U  ( I ) « (D ( I ) -VPROD (B ( 1 , 1) , NB , U , 1 , I -1 , XY ) ) /B ( I , 1 ) 

CHECK  VIOLATIONS;  FIND  ONE  MOST  NEGATIVE  LM 

XX-0,0 

K-0 

D05 J  *MP , NC 

XY  *VPROD  (C  ( 1 ,  J )  ,1,X,1,NX,XY)  -D  ( J ). 

IF(XY) 5,5,6 

IF(NP-NX)38,30,40 

STOP 

IF (NI ) 36,36,37 
D013I*1 , NI 

A  ( I ) * (VPROD (C ( 1 , 1 ) ,1,C(1,J) ,1,NX,XY)- 
L  VPROD ( A , 1 , B ( I , 1 ) ,NB,I-1,XY))/B(I,I) 

A (NP) -VPROD (C ( 1 , J) ,1,C(1,J),1,NX,XY) -VPROD { A , 1 , A , 1 , NI , XY) 
XU- (DfJ) -VPROD (A,1,U,1,NI ,XY) )/A (NP) 


'  '27  B  ( I  ,  J )  -  (VPROD  (0(1/1)  , ,C  ( 1 ,  J)  , 1 ,  NX , XY ) - 

1  VPROD <B( 1,1) , NB ,B  ( J  ,  1) ,  NB , J-l , XY) )/B(J, J) 

B (I, I) -SQRT (7PF0D(C(1,I) , 1 ,C ( 1 , I ) , 1 , NX , XY) - 
1  VPROD (B( 1,1) , MB , B ( 1 , l ) ,NB,IM,XY) ) 

26  U ( I ) ■ (D { I ) -VPROD (B ( 1 , 1) , NB, U , 1 , IM , XY) )  /B  ( I ,  I ) 

U(NI)-U(NI)/B(NI,NI) 

GO  TO  29 
17  0023J-1,NX 

I-UX+l-J 

23  X(I)n(X(I)+C(I,NC+l)-VPR0D(B(I,I+2) , NB , X ( 1+1 ) , 1 , J-l , XY) ) /B ( I , 1+1 ) 

IF (NI . EQ . 2 ) GOTO 30 
D012I-1,NI 
J-NI+l-I 

12  U (J) -(U(J) -VPROD (B(J+1,J) ,1,0<J41) ,1,1-1, XY) )/B(J,J) 

D039I-1 ,NC 

39  U  ( I ) -TN (LC ( I ) ) *U  (X ) 

30  IF(NP.GT.NC) G0T032 
D031J-NP,NC 
21  U(J)«0.0 

32  D033I-1,NC-1 
XX-U  (I) 

D03  4 J-l , NC 
IF (LC ( J) -I) 34,35,34 
.  v34  CONTINUE 
(£15  U(I)-U(J) 

U(J) -XX 
K-LC(I) 

LC  (I)  -LC (J) 

LC  ( J) -K 

33  CONTINUE 
RETURN 
END 

FUNCTION  VP, OD (A,L,B,M,N,  X) 

DIMENSION  A (1) ,B  (1) 

DOUBLE  PRECISION  XX 

XX-0.0 

IF (N) 1,1,2 

2  I-l-L 
J-l-M 
D03K-1 , N 
I-I+L 
J-J  +  M 

3  XX-VX  +  DBLE (A  a ) ) *DBLE (B (J) ) 

1  x-xx 

VPROD-XX 

RETURN 

END 

SUBROUTINE  MCDCMP ( A , L , N , X ) 

DIMENSION  A  ( L , 1 )  ,X(i) 

_  LE  (N)  3 , 3 , 4 
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It (XX-XU) 5,5,16 
.  16  K-J 

6$  XX-XU 

'*■'  A  (NP)  -SQRT  (A (NP)  } 

00191-1 ,NP 

19  B (NP , I ) -A (I ) 

5  CONTINUE 

IF  (K) 18,17,18 

C  IF  WE  HAVE  ONE,  COMPUTE  THE  LM ' S 

13  NI-NP 

NP-NP+1 

I-LC(NI) 

LC(NI) -LC(K) 

LC (K) -I 
XU-D(NI) 

D(NI)-D(K) 

D (K) -XU 
DO20I-1,NX 
XU-C ( I , NI ) 

C ( I , NI ) -C ( I , K) 

20  C ( I , K) -XU 
U (NI ) “XX 

29  XX-0.0 

K-0 

IF(NI.EQ.l) GOTO 10 
D021I -1 , NI-1 
J-NI-I 

V£r;  U  ( J)  ■  (U  ( J)  -VPROD  (B  ( J+l ,  J)  ,1,U(J+1)  ,  1, 1  ,XX)  )  /B  ( J,  J) 

IF (XX-U ( J) ) 22,21,21 
22  XX-U(J) 

J 

21  CONTINUE 

IF (K) 24,10,24 

C  IF  AN  LM>0  THEN  DELETE  THE  CONSTRAINT 

24  NP-NI 
NI -NI-1 
XX-D (NP) 

D (NP) -D (K) 

D(K) -XX 
I-LC (NP) 

LC (NP) -LC (K) 

LC  (K) -I 
D025I-1 , NX 
XX-C (I ,NP) 

C(I,NP)-C(I,K) 

25  C ( I , K ) "XX 

C  THEN  DO  THE  REDUCTION  FROM  K  ON 

D026I-K,‘NI 
IM-I-1 

IF (IM) 2b, 26, 8 
8  D027J -1 ,  IM 
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D01J-1,N 

JM-J-1 

A(J,J)-SQRT(A(J,J)-VPR0D(A(J,1) ,L,A(J,1) ,L,JMf2) ) 

IF  ( J .  EQ .  N) G0T01 

JP-J+1 

D02I-JP,N 

A ( X  f  J ) ■  ( A  ( I ,  J ) -VPROD ( A ( 1 , 1 ) ,L,A(J,1) ,L,JM,Z) )*X(J) 

CONTINUE 

RETURN 


Appendix  Ei  Comparison  of  Daughtar  Concantratlona 
Between  tha  Modified  Talvoglou  and 
laaat  Square  Mathoda 


Thia  appandlx  praaanta  In  tabular  form  tha  raaulta  of  tha  daughtar 
conoantrationa  calculated  uaing  both  mathoda  of  analysis.  Tha  tablaa  ara 
llatao  by  tha  catagory  of  runs  daacrlbad  in  Chaptar  V.  As  can  ba  saan  in 
tha  following  tablaa,  all  valuaa  axoapt  for  tha  conoantration  of  Ra-A  ara 
in  agraamant  batwaan  tha  two  mathoda.  Tha  values  of  Ra-A  in  all  caaas  ia 
off  by  a  factor  of  2,  which  indicataa  a  program  arror.  It  ia  raoommandad 
that  furthar  work  on  tha  program  ba  mads  to  datarmina  tha  aouroa  of  thia 
error  and  also  to  datarmina  sanaitivity  of  uaing  tha  laaat  square  method. 

TABLE  VZII 

Comparison  of  Valuaa  With  Blower  Only 


Taivoglou 

Weighted  Laaat 

Square 

RUN  1 

Ra-A 

Ra-B 

Ra-C  Ra-A  Ra-B 

Ra-C 

6.39 

4.52 

2.76 

1.79 

0.19 

0.89 

1.05 

2.54 

2.34 

2.22 

0.87 

2.94 

0.35 

0.41 

4.57 

3.90 

2.67 

2.98 
2.48 

1.99 

1.71 

1.02 

0.46 

1.59 

RUN  2 

4.34 

2.  B6 

3.67 

2.80 

2.05 

3.30 

2.09 

1.90 

2.05 

0.95 

0.71 

1.72 

1.58 

0,30 

2.63 

1.11 

0.64 

2.54 

1.24 

0.65 

2.14 

3.66 

3.64 


Taivoglou 


Weighted  Least  Square 


RUN  1 


Ra-A 

Ra-B 

Ra-C 

Ra-A 

Ra-B 

Ra-C 

3.27 

4.43 

2.94 

2.70 

3.87 

3.49 

3.52 

4.08 

2.86 

3.41 

4.01 

3.03 

1.93 

0.  11 

0.14 

1.30 

0.28 

0.11 

1.37 

0.04 

0.24 

1.09 

0.27 

0.17 

2.43 

0.40 

0.00 

1.38 

0.39 

0.02 

2.06 

0.11 

0.26 

1.02 

0.11 

0.38 

2.74 

0.34 

0.00 

1.26 

0.23 

0.10 

2.66 

0.24 

0.00 

1.32 

0.29 

0.02 

1 .  OB 

0.33 

4.01 

3.47 

3.23 

1.97 

6.83 

3.64 

2.15 

3.34 

3.62 

2.37 

RUN  2 

7.14 

4.39 

2.16 

3.70 

4.04 

2.83 

1.80 

0.72 

0.13 

0.97 

0.17 

0.17 

1.91 

0.22 

0.21 

0.91 

0.21 

0.33 

2.03 

0.35 

0.03 

1.43 

0.42 

0.03 

2.84 

0.  17 

0.00 

1.70 

0.31 

0.00 

2.38 

0. 13 

0.03 

1.33 

0.33 

0.02 

2.36 

0.23 

0.03 

1.21 

0.21 

0.19 

7.11 

2.42 

0.61 

3.77 

2.63 

0.93 

RUN  3 

7.32 

3.75 

3.95 

3.99 

3.03 

5.27 

2.94 

2.84 

3.68 

3.21 

2.19 

0.22 

0.02 

1.19 

0.24 

0.11 

2.21 

0.  14 

0.03 

1. 36 

0.24 

0.08 

2.32 

0.  16 

0.00 

1.17 

0.13 

0.34 

2.47 

0.18 

0.15 

1.53 

0.29 

0.65 

2.70 

0.24 

0.00 

1.78 

0.33 

0.00 

2.98 

0.28 

0.00 

1.50 

0.35 

0.03 

2.76 

0.27 

0.00 

1.69 

0.35 

0.02 

RUN  4 

Ml 

2.21 

0.09 

0.06 

1.36 

I 

0.08 

2.43 

0.15 

0.06 

1.66 

■3B1 

0.01 

0.99 

0.00 

0.51 

1.21 

0.24 

0.11 

2.47 

0.15 

0.02 

1.49 

0.23 

0.03 

2.36 

0.29 

0.00 

1.38 

0.37 

0.02 

2.13 

0.10 

0.03 

1.28 

0.21 

0.07 

TABLE  X 


Comparison  With  Half  of  Intaka  Saalad  on  tha  EAC 


Taivoglou 

Waightad  Laaat  Squara 

RUN  1 

Ra-B 

Ra-C  Ra»A 

Ra-B 

Ra-C 

4.64 

2.72 

2.18  3.18 

2.96 

2.11 

2.24 

0.32 

0.00  1.37 

0.45 

0.00 

2.55 

0.18 

0.00  1.52 

0.26 

0.01 

2.23 

0.18 

0.14  1.49 

0.37 

0.10 

2.48 

0.25 

0.00  1.36 

0.27 

0.10 

1.85 

0.10 

0.12  1.35 

0.33 

0.04 

2.42 

0.17 

0.00  1.46 

0.25 

0.02 

6.12 

2.09 

0.67  4.92 

3.58 

0.00 

5.75 

3.50 

2.02  3.94 

3.59 

2.11 

m 

RUN  2 

6.78 

2.99  5.19 

5.90 

4.72 

1.63 

0.86  2.13 

1.53 

1.13 

3.03 

0.16 

0.10  2.22 

0.47 

0.03 

4.04 

0.39 

0.00  2.08 

0.54 

0.00 

1.97 

0.09 

0.16  1.40 

0.32 

0.07 

1.89 

0.07 

0.09  1.17 

0.17 

0.08 

2.54 

0.21 

0.00  1.25 

0.21 

0.04 

5. 7e 

0.67 

0.06  3.33 

0.87 

0.21 

5.68 

2.99 

1.25  3.08 

2.87 

2.12 

RUN  3 

6.28 

3.95 

2.21  2.75 

3.43 

2.88 

2.25 

1.00 

0.53  1.62 

1.44 

0.45 

2.50 

0.03 

0.02  1.28 

0.23 

0.00 

2.60 

0.26 

0.00  1.21 

0.21 

0.09 

2.53 

0.24 

0.00  1.52 

0.36 

0.02 

2.33 

0.23 

0.00  1.66 

0.36 

0.00 

2.40 

0.19 

0.01  1.61 

0.35 

0.00 

5.87 

0.98 

0.00  3.42 

1.14 

0.14 

8.21 

3.41 

0.81  4.59 

3.71 

1. 18 
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